A micro introduction to Bayesian modeling with BRMS

Julin Maloof

Frequentist vs Bayesian Statistics

Frequentist Statistics

  • What is the probability of the data given a hypothesis
  • Theory relies on data that is not observed
  • Accept or reject null hypothesis
  • Does not require use of priors

Bayesian

  • What is the probability of the hypothesis given the data
  • Does not rely on unobserved data
  • Evaluate ratio of probabilities for different hypotheses
  • Requires specification of prior beliefs (not as bad as it sounds!)

Prior probabilities

  • Bayesian methods require that we specify prior probabilities on the model parameters
  • This might sound like we are cheating
  • Instead we use relatively uninformative priors
  • As long as we have set relatively weak (uninformative) priors, the data will overwhelm the priors
  • (Although if we do have prior information we can use informative priors)

Sampling

  • Bayesian model fitting requires sampling different parameter values
  • For each set of parameters test, calculate the probability of the model, given the data
    • Similar to what we discussed two weeks ago with the nls function
  • Difference is that for nls the goal was just to find the “BEST” set of parameters
  • For Bayesian we need to sample enough to define the probability “landscape”
  • Typically several thousand sets of parameter values are sampled

Posterior probabilities

  • One output of Bayesian model fitting is a set of posterior probabilities
  • These will give us the distribution of possible parameter values for each model parameter, given the data and the priors
  • We can do further tests/comparisons once we have these posterior probabilities

Diagnostics

  • It is critically important to evaluate whether the sampling worked well
  • Typically four sampling “chains” are run
  • Rhat compares variance between and within chains. Should be < 1.05
  • ESS (Effective sampling size) should be > 100 per chain

Let’s try it

library(tidyverse)
library(brms)
dat <- read_csv("../output/height_data_clean.csv")

Please create two objects

  • datsmall that only retains observations from “2023-01-27”
  • datsmallbcd that only retains BH, CC, and DPR observations from “2023-01-27”

Check the number of “cores” on our computer

  • Most computers have multiple “cores” or CPUs that can be used in parallel.
  • Let’s check what we have. We will use this on the next slide.
parallel::detectCores()
[1] 8

What are the default priors?

  • We will use what should now be a familiar formula to specify our model
  • We can ask what priors would be used by default for our formula and data
get_prior(height_cm ~ pop, 
          data = datsmallbcd)
                  prior     class   coef group resp dpar nlpar lb ub
                 (flat)         b                                   
                 (flat)         b  popCC                            
                 (flat)         b popDPR                            
 student_t(3, 3.5, 2.5) Intercept                                   
   student_t(3, 0, 2.5)     sigma                               0   
       source
      default
 (vectorized)
 (vectorized)
      default
      default
  • They “b” class are the coefficients on fixed effects.
  • They have a flat prior (all values from -Infinity to +Infinity are equally likely)
  • We can do better than this while still not biasing the model (next slide)
  • Student_t is like a fat-tailed normal distribution.
    • The intercept is centered at 3.5,
    • sigma (standard deviation) is centered at 0

Student vs normal

Fit a model with height being predicted by pop

m1 <- brm(height_cm ~ pop,
          data = datsmallbcd,
          prior = set_prior("normal(0,10)", class = "b"),
          sample_prior = TRUE,
          cores = 4) # set to 1 if you only have 1 core, or 2 if you have 2.
                       # don't go higher than 4, even if you have more than 4 cores.    
  • We use an uninformative prior for “b”, the coefficient for differences between pops
    • We are setting this prior to be normally distributed, centered on 0, with a standard deviation of 10
    • Note “b” refers to any fixed-effect coefficients. If we had multiple fixed effects we could give them all the same prior, or we could specify them individually.
  • By default brms uses 4 chains cores = 4 means that we want each chain to run in parallel in a different CPU in our computer.
  • By default, each chain has 1,000 warm-up samples and 1,000 real samples

Fit a model with height being predicted by pop

CORRECTED: cores=4 not threads=4

m1 <- brm(height_cm ~ pop,
          data = datsmallbcd,
          prior = set_prior("normal(0,10)", class = "b"),
          sample_prior = TRUE,
          cores = 4)

The output shows us the progress of each chain

Model summary

summary(m1)
 Family: gaussian 
  Links: mu = identity; sigma = identity 
Formula: height_cm ~ pop 
   Data: datsmallbcd (Number of observations: 155) 
  Draws: 4 chains, each with iter = 2000; warmup = 1000; thin = 1;
         total post-warmup draws = 4000

Population-Level Effects: 
          Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
Intercept     3.22      0.12     2.98     3.47 1.00     4278     2866
popCC         1.46      0.25     0.97     1.93 1.00     3895     3132
popDPR        1.91      0.31     1.30     2.51 1.00     4295     3206

Family Specific Parameters: 
      Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sigma     1.25      0.07     1.12     1.39 1.00     4186     2810

Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
and Tail_ESS are effective sample size measures, and Rhat is the potential
scale reduction factor on split chains (at convergence, Rhat = 1).
  • Estimate is the estimated height
    • Intercept is for the reference population (BH in this case)
    • popCC is the estimated difference between CC and BH
  • 95% confidence intervals show our confidence in the estimate.
    If they do not cross 0, then we are at least 95% confident that our differences are real
  • Be sure to check Rhat and the two ESS stats. What did we want these to be?

Model plots

plot(m1, N=3, ask=FALSE)

  • On the left we see the distribution of the posterior probability
    • Check for unusual pattern, bimodality, etc.
    • If the posterior is near the edge of the range of the prior, consider adjusting the prior
  • On the right we see the value for each parameter for each sample of each chain
    • Should look relatively uniform from left to right and chains should look similar

what are the other priors?

Our posterior plots had two coefficients that we did not specify priors for. What was used?

                  prior     class   coef group resp dpar nlpar lb ub
           normal(0,10)         b                                   
           normal(0,10)         b  popCC                            
           normal(0,10)         b popDPR                            
 student_t(3, 3.5, 2.5) Intercept                                   
   student_t(3, 0, 2.5)     sigma                               0   
       source
         user
 (vectorized)
 (vectorized)
      default
      default

Bayesian hypothesis testing

  • Is popCC really different from popBH?
  • Looking for
    • Post_prob to be near 1 (probability that hypothesis is true)
    • Evid.Ratio > 3 (>3 = “moderate evidence”, > 10 = “strong evidence”)
hypothesis(m1, hypothesis = "popCC > 0")
Hypothesis Tests for class b:
   Hypothesis Estimate Est.Error CI.Lower CI.Upper Evid.Ratio Post.Prob Star
1 (popCC) > 0     1.46      0.25     1.04     1.86        Inf         1    *
---
'CI': 90%-CI for one-sided and 95%-CI for two-sided hypotheses.
'*': For one-sided hypotheses, the posterior probability exceeds 95%;
for two-sided hypotheses, the value tested against lies outside the 95%-CI.
Posterior probabilities of point hypotheses assume equal prior probabilities.

Getting posterior probability for an estimate

Sometime I prefer to do this by “hand”

posterior <- m1 %>% as_draws_df()
dim(posterior)
[1] 4000   12
head(posterior)
# A draws_df: 6 iterations, 1 chains, and 9 variables
  b_Intercept b_popCC b_popDPR sigma prior_Intercept prior_b prior_sigma lprior
1         3.2     1.3      2.1   1.2             1.1    -6.0         5.3   -9.8
2         3.3     1.5      1.6   1.3             4.3    -6.1        29.5   -9.8
3         3.3     1.2      2.1   1.3             3.1     8.5         3.3   -9.8
4         3.2     1.6      2.3   1.2             5.3     9.3         2.2   -9.8
5         3.3     1.2      1.6   1.3            -5.7    16.0         3.8   -9.8
6         3.3     1.5      2.0   1.3             3.4    20.2         2.1   -9.8
# ... with 1 more variables
# ... hidden reserved variables {'.chain', '.iteration', '.draw'}

How many draws are less than 0?

posterior %>% 
  summarize(CC_greater_than_zero = sum(b_popCC > 0)/n())
# A tibble: 1 × 1
  CC_greater_than_zero
                 <dbl>
1                    1
  • All 4000 posterior samples estimate the the CC coefficient is > 0
  • This is very strong evidence the CC height is greater than BH height

Bayesian hypothesis testing

  • Is popCC the same height as popDPR?
  • Looking for
    • Post_prob to be near 1 (probability that hypothesis is true)
    • Evid.Ratio to > 3 (>3 = “moderate evidence”, > 10 = “strong evidence”)
hypothesis(m1, hypothesis = "popCC = popDPR")
Hypothesis Tests for class b:
            Hypothesis Estimate Est.Error CI.Lower CI.Upper Evid.Ratio
1 (popCC)-(popDPR) = 0    -0.45      0.36    -1.16     0.24      18.09
  Post.Prob Star
1      0.95     
---
'CI': 90%-CI for one-sided and 95%-CI for two-sided hypotheses.
'*': For one-sided hypotheses, the posterior probability exceeds 95%;
for two-sided hypotheses, the value tested against lies outside the 95%-CI.
Posterior probabilities of point hypotheses assume equal prior probabilities.

Fit a model with random effects:

  • Really pop and block should probably both be random, but we will keep pop fixed for now
m2 <- brm(height_cm ~ pop + (1|block),
          data = datsmallbcd,
          prior = set_prior("normal(0,10)", class = "b"),
          sample_prior = TRUE,
          threads = 4)

SAMPLING FOR MODEL 'anon_model' NOW (CHAIN 1).
Chain 1: 
Chain 1: Gradient evaluation took 0.000106 seconds
Chain 1: 1000 transitions using 10 leapfrog steps per transition would take 1.06 seconds.
Chain 1: Adjust your expectations accordingly!
Chain 1: 
Chain 1: 
Chain 1: Iteration:    1 / 2000 [  0%]  (Warmup)
Chain 1: Iteration:  200 / 2000 [ 10%]  (Warmup)
Chain 1: Iteration:  400 / 2000 [ 20%]  (Warmup)
Chain 1: Iteration:  600 / 2000 [ 30%]  (Warmup)
Chain 1: Iteration:  800 / 2000 [ 40%]  (Warmup)
Chain 1: Iteration: 1000 / 2000 [ 50%]  (Warmup)
Chain 1: Iteration: 1001 / 2000 [ 50%]  (Sampling)
Chain 1: Iteration: 1200 / 2000 [ 60%]  (Sampling)
Chain 1: Iteration: 1400 / 2000 [ 70%]  (Sampling)
Chain 1: Iteration: 1600 / 2000 [ 80%]  (Sampling)
Chain 1: Iteration: 1800 / 2000 [ 90%]  (Sampling)
Chain 1: Iteration: 2000 / 2000 [100%]  (Sampling)
Chain 1: 
Chain 1:  Elapsed Time: 0.187 seconds (Warm-up)
Chain 1:                0.13 seconds (Sampling)
Chain 1:                0.317 seconds (Total)
Chain 1: 

SAMPLING FOR MODEL 'anon_model' NOW (CHAIN 2).
Chain 2: 
Chain 2: Gradient evaluation took 1e-05 seconds
Chain 2: 1000 transitions using 10 leapfrog steps per transition would take 0.1 seconds.
Chain 2: Adjust your expectations accordingly!
Chain 2: 
Chain 2: 
Chain 2: Iteration:    1 / 2000 [  0%]  (Warmup)
Chain 2: Iteration:  200 / 2000 [ 10%]  (Warmup)
Chain 2: Iteration:  400 / 2000 [ 20%]  (Warmup)
Chain 2: Iteration:  600 / 2000 [ 30%]  (Warmup)
Chain 2: Iteration:  800 / 2000 [ 40%]  (Warmup)
Chain 2: Iteration: 1000 / 2000 [ 50%]  (Warmup)
Chain 2: Iteration: 1001 / 2000 [ 50%]  (Sampling)
Chain 2: Iteration: 1200 / 2000 [ 60%]  (Sampling)
Chain 2: Iteration: 1400 / 2000 [ 70%]  (Sampling)
Chain 2: Iteration: 1600 / 2000 [ 80%]  (Sampling)
Chain 2: Iteration: 1800 / 2000 [ 90%]  (Sampling)
Chain 2: Iteration: 2000 / 2000 [100%]  (Sampling)
Chain 2: 
Chain 2:  Elapsed Time: 0.171 seconds (Warm-up)
Chain 2:                0.15 seconds (Sampling)
Chain 2:                0.321 seconds (Total)
Chain 2: 

SAMPLING FOR MODEL 'anon_model' NOW (CHAIN 3).
Chain 3: 
Chain 3: Gradient evaluation took 1e-05 seconds
Chain 3: 1000 transitions using 10 leapfrog steps per transition would take 0.1 seconds.
Chain 3: Adjust your expectations accordingly!
Chain 3: 
Chain 3: 
Chain 3: Iteration:    1 / 2000 [  0%]  (Warmup)
Chain 3: Iteration:  200 / 2000 [ 10%]  (Warmup)
Chain 3: Iteration:  400 / 2000 [ 20%]  (Warmup)
Chain 3: Iteration:  600 / 2000 [ 30%]  (Warmup)
Chain 3: Iteration:  800 / 2000 [ 40%]  (Warmup)
Chain 3: Iteration: 1000 / 2000 [ 50%]  (Warmup)
Chain 3: Iteration: 1001 / 2000 [ 50%]  (Sampling)
Chain 3: Iteration: 1200 / 2000 [ 60%]  (Sampling)
Chain 3: Iteration: 1400 / 2000 [ 70%]  (Sampling)
Chain 3: Iteration: 1600 / 2000 [ 80%]  (Sampling)
Chain 3: Iteration: 1800 / 2000 [ 90%]  (Sampling)
Chain 3: Iteration: 2000 / 2000 [100%]  (Sampling)
Chain 3: 
Chain 3:  Elapsed Time: 0.183 seconds (Warm-up)
Chain 3:                0.157 seconds (Sampling)
Chain 3:                0.34 seconds (Total)
Chain 3: 

SAMPLING FOR MODEL 'anon_model' NOW (CHAIN 4).
Chain 4: 
Chain 4: Gradient evaluation took 1.2e-05 seconds
Chain 4: 1000 transitions using 10 leapfrog steps per transition would take 0.12 seconds.
Chain 4: Adjust your expectations accordingly!
Chain 4: 
Chain 4: 
Chain 4: Iteration:    1 / 2000 [  0%]  (Warmup)
Chain 4: Iteration:  200 / 2000 [ 10%]  (Warmup)
Chain 4: Iteration:  400 / 2000 [ 20%]  (Warmup)
Chain 4: Iteration:  600 / 2000 [ 30%]  (Warmup)
Chain 4: Iteration:  800 / 2000 [ 40%]  (Warmup)
Chain 4: Iteration: 1000 / 2000 [ 50%]  (Warmup)
Chain 4: Iteration: 1001 / 2000 [ 50%]  (Sampling)
Chain 4: Iteration: 1200 / 2000 [ 60%]  (Sampling)
Chain 4: Iteration: 1400 / 2000 [ 70%]  (Sampling)
Chain 4: Iteration: 1600 / 2000 [ 80%]  (Sampling)
Chain 4: Iteration: 1800 / 2000 [ 90%]  (Sampling)
Chain 4: Iteration: 2000 / 2000 [100%]  (Sampling)
Chain 4: 
Chain 4:  Elapsed Time: 0.188 seconds (Warm-up)
Chain 4:                0.145 seconds (Sampling)
Chain 4:                0.333 seconds (Total)
Chain 4: 

what priors were used for random effects?

prior_summary(m2)
                  prior     class      coef group resp dpar nlpar lb ub
           normal(0,10)         b                                      
           normal(0,10)         b     popCC                            
           normal(0,10)         b    popDPR                            
 student_t(3, 3.5, 2.5) Intercept                                      
   student_t(3, 0, 2.5)        sd                                  0   
   student_t(3, 0, 2.5)        sd           block                  0   
   student_t(3, 0, 2.5)        sd Intercept block                  0   
   student_t(3, 0, 2.5)     sigma                                  0   
       source
         user
 (vectorized)
 (vectorized)
      default
      default
 (vectorized)
 (vectorized)
      default

Check model sampling (do your own)

summary

 Family: gaussian 
  Links: mu = identity; sigma = identity 
Formula: height_cm ~ pop + (1 | block) 
   Data: datsmallbcd (Number of observations: 155) 
  Draws: 4 chains, each with iter = 2000; warmup = 1000; thin = 1;
         total post-warmup draws = 4000

Group-Level Effects: 
~block (Number of levels: 10) 
              Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sd(Intercept)     0.50      0.20     0.20     0.99 1.00     1119     1964

Population-Level Effects: 
          Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
Intercept     3.23      0.20     2.82     3.62 1.00     1623     1574
popCC         1.42      0.24     0.93     1.89 1.00     4132     2506
popDPR        1.86      0.30     1.28     2.43 1.00     3796     2604

Family Specific Parameters: 
      Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sigma     1.18      0.07     1.05     1.33 1.00     3864     3001

Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
and Tail_ESS are effective sample size measures, and Rhat is the potential
scale reduction factor on split chains (at convergence, Rhat = 1).

plots

Is it better than the model without the random effect?

m1 <- add_criterion(m1, "loo")
m2 <- add_criterion(m2, "loo")
loo_compare(m1, m2)
   elpd_diff se_diff
m2  0.0       0.0   
m1 -5.1       3.3   
  • Preferred model is listed first
  • elpd_diff is how much worse a model is relative to the preferred model
  • Want this to be greater (ideally 2X greater) than se_diff to favor the more complex model

Weibull in brms

Data set:

Use all dates for the Weibull modeling

dat <- read_csv("../output/height_data_clean.csv")

dat <- dat %>%
  mutate(day=as.numeric(survey_date-min(survey_date)+30),
         mf=factor(mf))

datbcd <- dat %>% filter(str_detect(pop, "BH|CC|DPR"))

Non-linear model specification in brms

Our formula specification is a bit more complex, so we use the function brmsformula to define the non-linear formula.

I am also going to rename a couple of the parameters to be more clear

f1 <- brmsformula(height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta),
                  Hmax + Hmin + k + delta ~ 1,
                  nl=TRUE)
  • The first line should look familiar (Weibull formula)
    • I am dividing k by 100 just to get it on a more reasonable scale
  • The second line is new, and it defines how each parameter should be fit with respect to the predictors in our data set.
  • The notation Hmax + Hmin + k + delta ~ 1 means that we want a single parameter value (Intercept) for each of the four parameters
  • nl=TRUE tells brms that this is a non-linear formula.

get the priors

Please review your notes and figure out how to get the possible priors for the formula on the previous slide

                prior class      coef group resp dpar nlpar lb ub       source
 student_t(3, 0, 4.2) sigma                                  0         default
               (flat)     b                           delta            default
               (flat)     b Intercept                 delta       (vectorized)
               (flat)     b                            Hmax            default
               (flat)     b Intercept                  Hmax       (vectorized)
               (flat)     b                            Hmin            default
               (flat)     b Intercept                  Hmin       (vectorized)
               (flat)     b                               k            default
               (flat)     b Intercept                     k       (vectorized)

Set the priors

What priors should we use?

prior1 <- c(set_prior("normal(5,3)", nlpar="Hmin"),
            set_prior("normal(60,15)",  nlpar="Hmax"),
            set_prior("gamma(2,2)", nlpar = "k"), # gamma constrains to be > 0
            set_prior("gamma(20, 3)", nlpar="delta")
)

Why gamma?

  • Want to constrain k and delta to be > 0.
  • Setting lower bounds causes problems
  • Gamma distribution is positive
tibble(x=seq(0,10,.1),
       gamma_2_2=dgamma(x, 2, 2),
       gamma_20_3=dgamma(x, 20, 3)) %>%
  pivot_longer(-x, values_to = "density") %>%
  ggplot(aes(x=x,y=density, color=name)) +
  geom_line()

Fit the model

fit1 <- brm(formula=f1, data=datbcd, prior=prior1, cores = 4)

Get the summary

summary(fit1)
 Family: gaussian 
  Links: mu = identity; sigma = identity 
Formula: height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta) 
         Hmax ~ 1
         Hmin ~ 1
         k ~ 1
         delta ~ 1
   Data: datbcd (Number of observations: 907) 
  Draws: 4 chains, each with iter = 2000; warmup = 1000; thin = 1;
         total post-warmup draws = 4000

Population-Level Effects: 
                Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
Hmax_Intercept     65.44      1.51    62.59    68.58 1.00     1877     1887
Hmin_Intercept      4.25      0.55     3.17     5.31 1.00     2909     3087
k_Intercept         0.67      0.01     0.64     0.69 1.00     1685     1932
delta_Intercept     6.12      0.49     5.22     7.13 1.01     1535     2176

Family Specific Parameters: 
      Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sigma    11.43      0.27    10.92    11.97 1.00     3342     2977

Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
and Tail_ESS are effective sample size measures, and Rhat is the potential
scale reduction factor on split chains (at convergence, Rhat = 1).

fit1 plots

plot(fit1, N = 3, ask=FALSE)

fit 1 predictions

fit1.prediction <- cbind(datbcd, pred1=predict(fit1)[,"Estimate"])

fit1.prediction %>% ggplot(aes(x=day)) +
  geom_line(aes(group=plantID, y=height_cm), alpha=.1) +
  geom_line(aes(y=pred1), color="blue") 

fit 1 plot predictions per pop

  • Please write code to generate the plot below (you can do this by adding one function call to the code on the previous slide)
  • In your own words, how did it do?

Allow a different Hmax for each pop {.nonincremental)}

  • The plants from different pops grow to different heights. Maybe we should model that.
  • With only 3 pops we will use a fixed effect. (Need 5+ for random effect)
  • See below for how we specify that Hmax should take on a different value for each pop
f2 <- brmsformula(height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta),
                  Hmin + k + delta ~ 1,
                  Hmax ~ pop,
                  nl=TRUE)

get the priors

Please review your notes and figure out how to get the possible priors for the formula on the previous slide

                prior class      coef group resp dpar nlpar lb ub       source
 student_t(3, 0, 4.2) sigma                                  0         default
               (flat)     b                           delta            default
               (flat)     b Intercept                 delta       (vectorized)
               (flat)     b                            Hmax            default
               (flat)     b Intercept                  Hmax       (vectorized)
               (flat)     b     popCC                  Hmax       (vectorized)
               (flat)     b    popDPR                  Hmax       (vectorized)
               (flat)     b                            Hmin            default
               (flat)     b Intercept                  Hmin       (vectorized)
               (flat)     b                               k            default
               (flat)     b Intercept                     k       (vectorized)

Set the priors

  • Do these seem reasonable to you?
prior2 <- c(set_prior("normal(5,3)", nlpar="Hmin"),
            set_prior("normal(60,15)", nlpar="Hmax"),
            set_prior("normal(0,20)", nlpar="Hmax", coef="popCC"),
            set_prior("normal(0,20)", nlpar="Hmax", coef="popDPR"),
            set_prior("gamma(2,2)", nlpar = "k"), # gamma constrains to be > 0
            set_prior("gamma(20, 3)", nlpar="delta")
)

Fit model 2

fit2 <- brm(formula=f2, data=datbcd, prior=prior2, cores = 4)

Fit 2 summary

Any Problems?

summary(fit2)
 Family: gaussian 
  Links: mu = identity; sigma = identity 
Formula: height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta) 
         Hmin ~ 1
         k ~ 1
         delta ~ 1
         Hmax ~ pop
   Data: datbcd (Number of observations: 907) 
  Draws: 4 chains, each with iter = 2000; warmup = 1000; thin = 1;
         total post-warmup draws = 4000

Population-Level Effects: 
                Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
Hmin_Intercept      4.01      0.53     2.93     5.02 1.00     3411     2577
k_Intercept         0.66      0.01     0.64     0.68 1.00     2617     2542
delta_Intercept     6.03      0.44     5.21     6.92 1.00     2644     2407
Hmax_Intercept     65.55      1.50    62.75    68.60 1.00     2829     2434
Hmax_popCC         10.61     11.70   -11.44    34.13 1.00     4229     2942
Hmax_popDPR        49.12     12.40    25.56    73.00 1.00     3972     2481

Family Specific Parameters: 
      Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sigma    11.28      0.27    10.77    11.82 1.00     3572     2405

Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
and Tail_ESS are effective sample size measures, and Rhat is the potential
scale reduction factor on split chains (at convergence, Rhat = 1).

Fit 2 plots

Any problems?

plot(fit2, N = 3, ask=FALSE)

Fit 2 plot the predictions

datbcd$pred2 <- predict(fit2)[,"Estimate"]

datbcd %>% ggplot(aes(x=day)) +
  geom_line(aes(group=plantID, y=height_cm, color=pop), alpha=.3) +
  geom_line(aes(y=pred2, color=pop)) +
  scale_color_brewer(type="qual", palette = "Set2", guide="none") +
  facet_wrap(~pop)

Fit 2 Extend the CC and DPR curves out for all days

  • It would be nice to get predictions for CC and DPR all the way out to day 200
  • To do this, we create a data frame containing each day and pop for which we would like a prediction
  • The function expand_grid creates a data frame with all combinations of the items it is given (days and pops in this case).
pred.df <- expand_grid(day=min(datbcd$day):max(datbcd$day),
                       pop=unique(datbcd$pop))
pred.df
# A tibble: 537 × 2
     day pop  
   <int> <chr>
 1    30 CC   
 2    30 BH   
 3    30 DPR  
 4    31 CC   
 5    31 BH   
 6    31 DPR  
 7    32 CC   
 8    32 BH   
 9    32 DPR  
10    33 CC   
# ℹ 527 more rows

Fit 2 Extend the CC and DPR curves out for all days (slide 2)

  • Now we use our prediction data frame to generate predictions
fit2.predictions <- cbind(pred.df, prediction=predict(fit2, newdata = pred.df)[,"Estimate"]) %>%
  full_join(datbcd, by=c("day", "pop"))

fit2.predictions
     day pop prediction plantID survey_date block row  col   mf rep height_cm
1     30  CC   4.179711    D14B  2023-01-27    D1   4    B    5  12       3.7
2     30  CC   4.179711    D16B  2023-01-27    D1   6    B    3  10       2.3
3     30  CC   4.179711   D112A  2023-01-27    D1  12    A    7   4       4.1
4     30  CC   4.179711   D228B  2023-01-27    D2  28    B    2  13       4.1
5     30  CC   4.179711   D235A  2023-01-27    D2  35    A    1  11       5.4
6     30  CC   4.179711   D239A  2023-01-27    D2  39    A    5   7       3.6
7     30  CC   4.179711   D239B  2023-01-27    D2  39    B    4   5       5.1
8     30  CC   4.179711    D14C  2023-01-27    D1   4    C    4   9       0.9
9     30  CC   4.179711   F112A  2023-01-27    F1  12    A    6   3       6.5
10    30  CC   4.179711   F120B  2023-01-27    F1  20    B    2   4       6.0
11    30  CC   4.179711   F231A  2023-01-27    F2  31    A    7   8       4.9
12    30  CC   4.179711   F117D  2023-01-27    F1  17    D    3  13       4.7
13    30  CC   4.179711   F224C  2023-01-27    F2  24    C    1   6       3.4
14    30  CC   4.179711   F234D  2023-01-27    F2  34    D    3  12       5.6
15    30  CC   4.179711   H111B  2023-01-27    H1  11    B    7   2       6.8
16    30  CC   4.179711   H115B  2023-01-27    H1  15    B    3   3       5.5
17    30  CC   4.179711   H116A  2023-01-27    H1  16    A    4  10       6.3
18    30  CC   4.179711   H116B  2023-01-27    H1  16    B    6   6       4.4
19    30  CC   4.179711   H221A  2023-01-27    H2  21    A    3   5       4.8
20    30  CC   4.179711   H234B  2023-01-27    H2  34    B    4   7       3.4
21    30  CC   4.179711    H19D  2023-01-27    H1   9    D    2   2       4.8
22    30  CC   4.179711   H236C  2023-01-27    H2  36    C    2   1       5.6
23    30  CC   4.179711    J19D  2023-01-27    J1   9    D    2   6       3.8
24    30  CC   4.179711   J116C  2023-01-27    J1  16    C    6   4       5.6
25    30  CC   4.179711   J236C  2023-01-27    J2  36    C    5   5       4.2
26    30  CC   4.179711   J239C  2023-01-27    J2  39    C    2   3       6.9
27    30  CC   4.179711   J110A  2023-01-27    J1  10    A    3   2       6.3
28    30  CC   4.179711   J112B  2023-01-27    J1  12    B    4   4       4.1
29    30  CC   4.179711   J238B  2023-01-27    J2  38    B    3   9       4.6
30    30  CC   4.179711    L18A  2023-01-27    L1   8    A    4   6       3.4
31    30  CC   4.179711   L232A  2023-01-27    L2  32    A    5   3       5.5
32    30  CC   4.179711   L114D  2023-01-27    L1  14    D    1   1       3.9
33    30  CC   4.179711   L117C  2023-01-27    L1  17    C    7   3       4.2
34    30  BH   3.803480    D15B  2023-01-27    D1   5    B    5  24       3.9
35    30  BH   3.803480    D18B  2023-01-27    D1   8    B    6   7       2.2
36    30  BH   3.803480   D111B  2023-01-27    D1  11    B    7   4       3.7
37    30  BH   3.803480   D119B  2023-01-27    D1  19    B    2   8       3.7
38    30  BH   3.803480   D120A  2023-01-27    D1  20    A    4   2       3.3
39    30  BH   3.803480   D226B  2023-01-27    D2  26    B    5  15       2.7
40    30  BH   3.803480   D230B  2023-01-27    D2  30    B    4   3       3.0
41    30  BH   3.803480   D233B  2023-01-27    D2  33    B    6   8       2.2
42    30  BH   3.803480   D235B  2023-01-27    D2  35    B    5  25       2.6
43    30  BH   3.803480    D16C  2023-01-27    D1   6    C    4  12       3.1
44    30  BH   3.803480    D18D  2023-01-27    D1   8    D    7  17       1.6
45    30  BH   3.803480    D19C  2023-01-27    D1   9    C    1   8       2.1
46    30  BH   3.803480   D117C  2023-01-27    D1  17    C    5  14       3.2
47    30  BH   3.803480   D118D  2023-01-27    D1  18    D    3   5       2.3
48    30  BH   3.803480   D119D  2023-01-27    D1  19    D    5   9       1.4
49    30  BH   3.803480   D222C  2023-01-27    D2  22    C    7  18       1.4
50    30  BH   3.803480   D227C  2023-01-27    D2  27    C    7   5       3.1
51    30  BH   3.803480   D229C  2023-01-27    D2  29    C    4  13       3.3
52    30  BH   3.803480   D229D  2023-01-27    D2  29    D    2   9       2.9
53    30  BH   3.803480   D230D  2023-01-27    D2  30    D    1  10       2.3
54    30  BH   3.803480   D235D  2023-01-27    D2  35    D    3   6       3.7
55    30  BH   3.803480   F110A  2023-01-27    F1  10    A    7  19       2.9
56    30  BH   3.803480   F115A  2023-01-27    F1  15    A    6   9       3.6
57    30  BH   3.803480   F117B  2023-01-27    F1  17    B    2  10       3.4
58    30  BH   3.803480   F227A  2023-01-27    F2  27    A    7   7       3.4
59    30  BH   3.803480   F229B  2023-01-27    F2  29    B    4   5       4.8
60    30  BH   3.803480   F240A  2023-01-27    F2  40    A    3   8       3.6
61    30  BH   3.803480   F240B  2023-01-27    F2  40    B    5  12       6.5
62    30  BH   3.803480    F13C  2023-01-27    F1   3    C    5  16       2.9
63    30  BH   3.803480    F15C  2023-01-27    F1   5    C    5  26       2.9
64    30  BH   3.803480    F17D  2023-01-27    F1   7    D    7   6       3.5
65    30  BH   3.803480   F113C  2023-01-27    F1  13    C    3   7       2.8
66    30  BH   3.803480   F115C  2023-01-27    F1  15    C    4   4       3.4
67    30  BH   3.803480   F120C  2023-01-27    F1  20    C    5  11       5.1
68    30  BH   3.803480   F120D  2023-01-27    F1  20    D    5   1       2.9
69    30  BH   3.803480   F225C  2023-01-27    F2  25    C    7  20       2.4
70    30  BH   3.803480   F231C  2023-01-27    F2  31    C    6  10       3.0
71    30  BH   3.803480   F233C  2023-01-27    F2  33    C    5  17       3.0
72    30  BH   3.803480   F233D  2023-01-27    F2  33    D    2  11       3.6
73    30  BH   3.803480   F235C  2023-01-27    F2  35    C    5  27       3.2
74    30  BH   3.803480   F235D  2023-01-27    F2  35    D    2   1       3.7
75    30  BH   3.803480    H14B  2023-01-27    H1   4    B    2  12       4.5
76    30  BH   3.803480    H18B  2023-01-27    H1   8    B    3   9       3.2
77    30  BH   3.803480   H114A  2023-01-27    H1  14    A    5   3       6.9
78    30  BH   3.803480   H117A  2023-01-27    H1  17    A    4   6       2.3
79    30  BH   3.803480   H118B  2023-01-27    H1  18    B    5  18       3.5
80    30  BH   3.803480   H224B  2023-01-27    H2  24    B    5  29       4.5
81    30  BH   3.803480   H227A  2023-01-27    H2  27    A    5   4       4.5
82    30  BH   3.803480   H229B  2023-01-27    H2  29    B    5  19       3.2
83    30  BH   3.803480   H233A  2023-01-27    H2  33    A    3  10       2.5
84    30  BH   3.803480   H233B  2023-01-27    H2  33    B    4   7       2.9
85    30  BH   3.803480   H239A  2023-01-27    H2  39    A    7   9       3.0
86    30  BH   3.803480   H239B  2023-01-27    H2  39    B    6   2       3.5
87    30  BH   3.803480    H15C  2023-01-27    H1   5    C    5  28       4.6
88    30  BH   3.803480   H110C  2023-01-27    H1  10    C    6  11       4.1
89    30  BH   3.803480   H115D  2023-01-27    H1  15    D    5  13       4.1
90    30  BH   3.803480   H117C  2023-01-27    H1  17    C    2   2       3.5
91    30  BH   3.803480   H119C  2023-01-27    H1  19    C    7  21       4.5
92    30  BH   3.803480   H222D  2023-01-27    H2  22    D    2   3       2.8
93    30  BH   3.803480   H223C  2023-01-27    H2  23    C    7  22       2.8
94    30  BH   3.803480   H227D  2023-01-27    H2  27    D    6  12       2.8
95    30  BH   3.803480   H235C  2023-01-27    H2  35    C    2  13       2.5
96    30  BH   3.803480    J19C  2023-01-27    J1   9    C    5  30       4.4
97    30  BH   3.803480   J110C  2023-01-27    J1  10    C    3   1       1.7
98    30  BH   3.803480   J112C  2023-01-27    J1  12    C    7  23       2.0
99    30  BH   3.803480   J114D  2023-01-27    J1  14    D    4   8       3.3
100   30  BH   3.803480   J116D  2023-01-27    J1  16    D    5   5       2.7
101   30  BH   3.803480   J119C  2023-01-27    J1  19    C    5  20       4.3
102   30  BH   3.803480   J223D  2023-01-27    J2  23    D    4   9       4.3
103   30  BH   3.803480   J225D  2023-01-27    J2  25    D    2   5       1.9
104   30  BH   3.803480   J230C  2023-01-27    J2  30    C    4  14       4.4
105   30  BH   3.803480   J231C  2023-01-27    J2  31    C    5   6       6.1
106   30  BH   3.803480   J231D  2023-01-27    J2  31    D    7  14       2.7
107   30  BH   3.803480   J232D  2023-01-27    J2  32    D    3  12       3.5
108   30  BH   3.803480   J235C  2023-01-27    J2  35    C    7  11       1.5
109   30  BH   3.803480   J236D  2023-01-27    J2  36    D    5  21       3.0
110   30  BH   3.803480    J14B  2023-01-27    J1   4    B    2   4       2.3
111   30  BH   3.803480   J119B  2023-01-27    J1  19    B    3  11       3.8
112   30  BH   3.803480   J120B  2023-01-27    J1  20    B    7  10       2.2
113   30  BH   3.803480   J241B  2023-01-27    J2  41    B    6   4       5.5
114   30  BH   3.803480   J233A  2023-01-27    J2  33    A    3   2       3.0
115   30  BH   3.803480   J237B  2023-01-27    J2  37    B    1   7       3.0
116   30  BH   3.803480    L15B  2023-01-27    L1   5    B    2   6       2.8
117   30  BH   3.803480    L18B  2023-01-27    L1   8    B    7   2       3.3
118   30  BH   3.803480   L110A  2023-01-27    L1  10    A    2  14       1.4
119   30  BH   3.803480   L111B  2023-01-27    L1  11    B    3   3       2.7
120   30  BH   3.803480   L115A  2023-01-27    L1  15    A    5   7       3.0
121   30  BH   3.803480   L120A  2023-01-27    L1  20    A    7  15       2.3
122   30  BH   3.803480    L14D  2023-01-27    L1   4    D    4  10       2.1
123   30  BH   3.803480   L113C  2023-01-27    L1  13    C    3  13       2.4
124   30  BH   3.803480   L113D  2023-01-27    L1  13    D    5  22       3.7
125   30  BH   3.803480   L116C  2023-01-27    L1  16    C    6   5       4.4
126   30  BH   3.803480   L118D  2023-01-27    L1  18    D    7  12       3.7
127   30  BH   3.803480   L228C  2023-01-27    L2  28    C    6   6       3.1
128   30  BH   3.803480   L228D  2023-01-27    L2  28    D    2   7       2.9
129   30  BH   3.803480   L229D  2023-01-27    L2  29    D    3  14       1.2
130   30  BH   3.803480   L230C  2023-01-27    L2  30    C    4  11       4.5
131   30  BH   3.803480   L234C  2023-01-27    L2  34    C    5   8       3.6
132   30  BH   3.803480   L237D  2023-01-27    L2  37    D    7  16       2.6
133   30  BH   3.803480   L240D  2023-01-27    L2  40    D    7  13       3.8
134   30  BH   3.803480   L241D  2023-01-27    L2  41    D    4   1       3.3
135   30  BH   3.803480   L242C  2023-01-27    L2  42    C    7   3       2.0
136   30 DPR   4.123882   D117A  2023-01-27    D1  17    A    6   4       4.1
137   30 DPR   4.123882   D225A  2023-01-27    D2  25    A    2   7       8.5
138   30 DPR   4.123882   D227B  2023-01-27    D2  27    B    4   9       7.0
139   30 DPR   4.123882    F19A  2023-01-27    F1   9    A    2   5       6.5
140   30 DPR   4.123882   F232A  2023-01-27    F2  32    A    6   2       2.5
141   30 DPR   4.123882   F241A  2023-01-27    F2  41    A    2   4       5.5
142   30 DPR   4.123882    F16C  2023-01-27    F1   6    C    6   9       2.5
143   30 DPR   4.123882   F240C  2023-01-27    F2  40    C    7   4       5.4
144   30 DPR   4.123882   H117B  2023-01-27    H1  17    B    5   2       6.4
145   30 DPR   4.123882   H226A  2023-01-27    H2  26    A    4  13       3.5
146   30 DPR   4.123882    H17C  2023-01-27    H1   7    C    2  10       7.8
147   30 DPR   4.123882   H225C  2023-01-27    H2  25    C    2   6       7.5
148   30 DPR   4.123882   J121D  2023-01-27    J1  21    D    4   8       3.9
149   30 DPR   4.123882    J16B  2023-01-27    J1   6    B    6  10       4.1
150   30 DPR   4.123882   J232A  2023-01-27    J2  32    A    2  12       7.5
151   30 DPR   4.123882   J232B  2023-01-27    J2  32    B    6  11       3.2
152   30 DPR   4.123882   L119B  2023-01-27    L1  19    B    4  11       4.3
153   30 DPR   4.123882   L237B  2023-01-27    L2  37    B    6   3       2.3
154   30 DPR   4.123882   L116D  2023-01-27    L1  16    D    6  13       4.5
155   30 DPR   4.123882   L233C  2023-01-27    L2  33    C    2   2       5.4
156   31  CC   4.134878    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
157   31  BH   4.285610    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
158   31 DPR   3.998382    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
159   32  CC   4.039877    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
160   32  BH   4.023100    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
161   32 DPR   3.894736    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
162   33  CC   4.055874    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
163   33  BH   3.882614    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
164   33 DPR   3.847568    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
165   34  CC   4.475893    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
166   34  BH   4.031812    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
167   34 DPR   4.133336    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
168   35  CC   3.947568    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
169   35  BH   3.979473    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
170   35 DPR   3.855128    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
171   36  CC   3.860387    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
172   36  BH   4.384213    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
173   36 DPR   3.995199    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
174   37  CC   3.894625    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
175   37  BH   4.429385    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
176   37 DPR   4.176377    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
177   38  CC   3.967851    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
178   38  BH   4.213488    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
179   38 DPR   4.189134    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
180   39  CC   3.981096    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
181   39  BH   3.911712    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
182   39 DPR   4.291499    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
183   40  CC   4.143798    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
184   40  BH   4.216873    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
185   40 DPR   3.956775    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
186   41  CC   3.951514    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
187   41  BH   4.072672    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
188   41 DPR   3.698380    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
189   42  CC   4.515231    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
190   42  BH   3.892124    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
191   42 DPR   3.833676    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
192   43  CC   3.916465    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
193   43  BH   4.103789    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
194   43 DPR   4.262297    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
195   44  CC   4.154131    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
196   44  BH   3.830371    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
197   44 DPR   3.992415    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
198   45  CC   4.162432    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
199   45  BH   4.007781    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
200   45 DPR   4.315733    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
201   46  CC   4.210873    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
202   46  BH   4.270516    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
203   46 DPR   4.169308    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
204   47  CC   4.311016    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
205   47  BH   4.164411    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
206   47 DPR   4.055664    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
207   48  CC   4.216497    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
208   48  BH   3.831159    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
209   48 DPR   4.308045    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
210   49  CC   4.020645    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
211   49  BH   4.167535    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
212   49 DPR   4.056207    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
213   50  CC   4.081845    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
214   50  BH   3.889019    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
215   50 DPR   3.896542    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
216   51  CC   4.019384    D14B  2023-02-17    D1   4    B    5  12       3.4
217   51  CC   4.019384   D112A  2023-02-17    D1  12    A    7   4       4.3
218   51  CC   4.019384   D228B  2023-02-17    D2  28    B    2  13       4.0
219   51  CC   4.019384   D235A  2023-02-17    D2  35    A    1  11       5.5
220   51  CC   4.019384   D239A  2023-02-17    D2  39    A    5   7       3.8
221   51  CC   4.019384   D239B  2023-02-17    D2  39    B    4   5       5.2
222   51  CC   4.019384   F112A  2023-02-17    F1  12    A    6   3       6.1
223   51  CC   4.019384   F120B  2023-02-17    F1  20    B    2   4       5.6
224   51  CC   4.019384   F231A  2023-02-17    F2  31    A    7   8       4.1
225   51  CC   4.019384   F117D  2023-02-17    F1  17    D    3  13       4.6
226   51  CC   4.019384   F234D  2023-02-17    F2  34    D    3  12       5.5
227   51  CC   4.019384   H111B  2023-02-17    H1  11    B    7   2       7.6
228   51  CC   4.019384   H115B  2023-02-17    H1  15    B    3   3       6.1
229   51  CC   4.019384   H116A  2023-02-17    H1  16    A    4  10       6.1
230   51  CC   4.019384   H234B  2023-02-17    H2  34    B    4   7       2.9
231   51  CC   4.019384    H19D  2023-02-17    H1   9    D    2   2       4.0
232   51  CC   4.019384   H236C  2023-02-17    H2  36    C    2   1       6.2
233   51  CC   4.019384   J110A  2023-02-17    J1  10    A    3   2       6.5
234   51  CC   4.019384   J112B  2023-02-17    J1  12    B    4   4       4.5
235   51  CC   4.019384   J228A  2023-02-17    J2  28    A    5   6       3.6
236   51  CC   4.019384   J238B  2023-02-17    J2  38    B    3   9       3.5
237   51  CC   4.019384    J19D  2023-02-17    J1   9    D    2   6       3.1
238   51  CC   4.019384   J236C  2023-02-17    J2  36    C    5   5       4.9
239   51  CC   4.019384   J239C  2023-02-17    J2  39    C    2   3       7.7
240   51  CC   4.019384    L18A  2023-02-17    L1   8    A    4   6       4.1
241   51  CC   4.019384   L232A  2023-02-17    L2  32    A    5   3       5.9
242   51  CC   4.019384   L117C  2023-02-17    L1  17    C    7   3       4.9
243   51  BH   3.968905    D15B  2023-02-17    D1   5    B    5  24       4.3
244   51  BH   3.968905    D18B  2023-02-17    D1   8    B    6   7       2.9
245   51  BH   3.968905   D111B  2023-02-17    D1  11    B    7   4       3.5
246   51  BH   3.968905   D119B  2023-02-17    D1  19    B    2   8       3.3
247   51  BH   3.968905   D120A  2023-02-17    D1  20    A    4   2       3.9
248   51  BH   3.968905   D226B  2023-02-17    D2  26    B    5  15       3.0
249   51  BH   3.968905   D230B  2023-02-17    D2  30    B    4   3       3.9
250   51  BH   3.968905   D233B  2023-02-17    D2  33    B    6   8       2.2
251   51  BH   3.968905   D235B  2023-02-17    D2  35    B    5  25       1.9
252   51  BH   3.968905    D16C  2023-02-17    D1   6    C    4  12       4.7
253   51  BH   3.968905    D18D  2023-02-17    D1   8    D    7  17       1.5
254   51  BH   3.968905    D19C  2023-02-17    D1   9    C    1   8      24.0
255   51  BH   3.968905   D117C  2023-02-17    D1  17    C    5  14       2.9
256   51  BH   3.968905   D118D  2023-02-17    D1  18    D    3   5       2.2
257   51  BH   3.968905   D119D  2023-02-17    D1  19    D    5   9       0.9
258   51  BH   3.968905   D227C  2023-02-17    D2  27    C    7   5       3.4
259   51  BH   3.968905   D229C  2023-02-17    D2  29    C    4  13       3.4
260   51  BH   3.968905   D229D  2023-02-17    D2  29    D    2   9       3.3
261   51  BH   3.968905   D230D  2023-02-17    D2  30    D    1  10       2.1
262   51  BH   3.968905   D235D  2023-02-17    D2  35    D    3   6       4.0
263   51  BH   3.968905   F115A  2023-02-17    F1  15    A    6   9       4.1
264   51  BH   3.968905   F117B  2023-02-17    F1  17    B    2  10       3.0
265   51  BH   3.968905   F229B  2023-02-17    F2  29    B    4   5       5.1
266   51  BH   3.968905   F240B  2023-02-17    F2  40    B    5  12       4.2
267   51  BH   3.968905    F15C  2023-02-17    F1   5    C    5  26       2.9
268   51  BH   3.968905    F17D  2023-02-17    F1   7    D    7   6       3.1
269   51  BH   3.968905   F113C  2023-02-17    F1  13    C    3   7       2.9
270   51  BH   3.968905   F115C  2023-02-17    F1  15    C    4   4       3.5
271   51  BH   3.968905   F120D  2023-02-17    F1  20    D    5   1       1.7
272   51  BH   3.968905   F225C  2023-02-17    F2  25    C    7  20       2.2
273   51  BH   3.968905   F233C  2023-02-17    F2  33    C    5  17       3.1
274   51  BH   3.968905   F233D  2023-02-17    F2  33    D    2  11       3.7
275   51  BH   3.968905   F235C  2023-02-17    F2  35    C    5  27       3.5
276   51  BH   3.968905   F235D  2023-02-17    F2  35    D    2   1       3.6
277   51  BH   3.968905    H14B  2023-02-17    H1   4    B    2  12       5.2
278   51  BH   3.968905    H18B  2023-02-17    H1   8    B    3   9       3.3
279   51  BH   3.968905   H114A  2023-02-17    H1  14    A    5   3       7.4
280   51  BH   3.968905   H117A  2023-02-17    H1  17    A    4   6       3.5
281   51  BH   3.968905   H118B  2023-02-17    H1  18    B    5  18       3.2
282   51  BH   3.968905   H224B  2023-02-17    H2  24    B    5  29       4.6
283   51  BH   3.968905   H227A  2023-02-17    H2  27    A    5   4       3.9
284   51  BH   3.968905   H229B  2023-02-17    H2  29    B    5  19       3.1
285   51  BH   3.968905   H233B  2023-02-17    H2  33    B    4   7       2.5
286   51  BH   3.968905   H239B  2023-02-17    H2  39    B    6   2       3.9
287   51  BH   3.968905    H15C  2023-02-17    H1   5    C    5  28       4.7
288   51  BH   3.968905   H110C  2023-02-17    H1  10    C    6  11       4.0
289   51  BH   3.968905   H115D  2023-02-17    H1  15    D    5  13       4.5
290   51  BH   3.968905   H119C  2023-02-17    H1  19    C    7  21       4.1
291   51  BH   3.968905   H222D  2023-02-17    H2  22    D    2   3       2.4
292   51  BH   3.968905   H223C  2023-02-17    H2  23    C    7  22       3.2
293   51  BH   3.968905   H227D  2023-02-17    H2  27    D    6  12       3.8
294   51  BH   3.968905   H235C  2023-02-17    H2  35    C    2  13       2.2
295   51  BH   3.968905    J14B  2023-02-17    J1   4    B    2   4       2.6
296   51  BH   3.968905   J119B  2023-02-17    J1  19    B    3  11       4.2
297   51  BH   3.968905   J233A  2023-02-17    J2  33    A    3   2       3.1
298   51  BH   3.968905   J237B  2023-02-17    J2  37    B    1   7       2.7
299   51  BH   3.968905   J241B  2023-02-17    J2  41    B    6   4       4.5
300   51  BH   3.968905    J19C  2023-02-17    J1   9    C    5  30       5.5
301   51  BH   3.968905   J110C  2023-02-17    J1  10    C    3   1       1.9
302   51  BH   3.968905   J114D  2023-02-17    J1  14    D    4   8       3.4
303   51  BH   3.968905   J119C  2023-02-17    J1  19    C    5  20       5.3
304   51  BH   3.968905   J223D  2023-02-17    J2  23    D    4   9       1.7
305   51  BH   3.968905   J225D  2023-02-17    J2  25    D    2   5       1.9
306   51  BH   3.968905   J230C  2023-02-17    J2  30    C    4  14       5.3
307   51  BH   3.968905   J231C  2023-02-17    J2  31    C    5   6       7.1
308   51  BH   3.968905   J231D  2023-02-17    J2  31    D    7  14       4.8
309   51  BH   3.968905   J232D  2023-02-17    J2  32    D    3  12       3.4
310   51  BH   3.968905   J236D  2023-02-17    J2  36    D    5  21      25.0
311   51  BH   3.968905    L15B  2023-02-17    L1   5    B    2   6       2.9
312   51  BH   3.968905    L18B  2023-02-17    L1   8    B    7   2       2.7
313   51  BH   3.968905   L110A  2023-02-17    L1  10    A    2  14       0.6
314   51  BH   3.968905   L111B  2023-02-17    L1  11    B    3   3       2.1
315   51  BH   3.968905   L120A  2023-02-17    L1  20    A    7  15       2.3
316   51  BH   3.968905    L14D  2023-02-17    L1   4    D    4  10       2.2
317   51  BH   3.968905   L113C  2023-02-17    L1  13    C    3  13       2.2
318   51  BH   3.968905   L113D  2023-02-17    L1  13    D    5  22       4.7
319   51  BH   3.968905   L116C  2023-02-17    L1  16    C    6   5       4.4
320   51  BH   3.968905   L228C  2023-02-17    L2  28    C    6   6       2.4
321   51  BH   3.968905   L228D  2023-02-17    L2  28    D    2   7       2.8
322   51  BH   3.968905   L229D  2023-02-17    L2  29    D    3  14       1.9
323   51  BH   3.968905   L230C  2023-02-17    L2  30    C    4  11       4.4
324   51  BH   3.968905   L234C  2023-02-17    L2  34    C    5   8       4.1
325   51  BH   3.968905   L237D  2023-02-17    L2  37    D    7  16       3.6
326   51  BH   3.968905   L240D  2023-02-17    L2  40    D    7  13       4.1
327   51  BH   3.968905   L241D  2023-02-17    L2  41    D    4   1       3.5
328   51 DPR   4.292258   D225A  2023-02-17    D2  25    A    2   7       8.9
329   51 DPR   4.292258   D227B  2023-02-17    D2  27    B    4   9       7.2
330   51 DPR   4.292258    F19A  2023-02-17    F1   9    A    2   5       6.0
331   51 DPR   4.292258   F232A  2023-02-17    F2  32    A    6   2       2.0
332   51 DPR   4.292258    F16C  2023-02-17    F1   6    C    6   9       2.6
333   51 DPR   4.292258   F240C  2023-02-17    F2  40    C    7   4       6.3
334   51 DPR   4.292258   H117B  2023-02-17    H1  17    B    5   2       8.1
335   51 DPR   4.292258    H17C  2023-02-17    H1   7    C    2  10       8.2
336   51 DPR   4.292258   H225C  2023-02-17    H2  25    C    2   6       7.8
337   51 DPR   4.292258    J16B  2023-02-17    J1   6    B    6  10       3.6
338   51 DPR   4.292258   J232B  2023-02-17    J2  32    B    6  11       2.1
339   51 DPR   4.292258   L119B  2023-02-17    L1  19    B    4  11       3.7
340   51 DPR   4.292258   L237B  2023-02-17    L2  37    B    6   3       4.0
341   51 DPR   4.292258   L116D  2023-02-17    L1  16    D    6  13       4.0
342   51 DPR   4.292258   L233C  2023-02-17    L2  33    C    2   2       5.7
343   52  CC   3.817100    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
344   52  BH   4.178649    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
345   52 DPR   4.109877    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
346   53  CC   3.751717    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
347   53  BH   4.025464    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
348   53 DPR   4.241984    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
349   54  CC   3.951848    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
350   54  BH   3.807669    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
351   54 DPR   4.376220    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
352   55  CC   4.275465    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
353   55  BH   4.158485    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
354   55 DPR   4.178879    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
355   56  CC   3.911960    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
356   56  BH   4.132348    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
357   56 DPR   4.331887    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
358   57  CC   4.097711    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
359   57  BH   4.413554    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
360   57 DPR   4.356327    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
361   58  CC   4.220583    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
362   58  BH   4.210526    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
363   58 DPR   4.281331    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
364   59  CC   4.406154    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
365   59  BH   4.420057    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
366   59 DPR   4.484808    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
367   60  CC   3.950863    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
368   60  BH   4.043791    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
369   60 DPR   4.685506    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
370   61  CC   4.146779    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
371   61  BH   4.322911    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
372   61 DPR   4.273247    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
373   62  CC   4.155638    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
374   62  BH   4.191693    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
375   62 DPR   4.819814    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
376   63  CC   4.488627    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
377   63  BH   4.461917    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
378   63 DPR   4.890417    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
379   64  CC   4.244708    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
380   64  BH   4.314203    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
381   64 DPR   4.358447    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
382   65  CC   4.594647   D112A  2023-03-03    D1  12    A    7   4       5.7
383   65  CC   4.594647   D235A  2023-03-03    D2  35    A    1  11       5.9
384   65  CC   4.594647   D239A  2023-03-03    D2  39    A    5   7       3.9
385   65  CC   4.594647   D239B  2023-03-03    D2  39    B    4   5       5.8
386   65  CC   4.594647   F231A  2023-03-03    F2  31    A    7   8       5.4
387   65  CC   4.594647   F117D  2023-03-03    F1  17    D    3  13       5.7
388   65  CC   4.594647   H111B  2023-03-03    H1  11    B    7   2       8.6
389   65  CC   4.594647   H115B  2023-03-03    H1  15    B    3   3       6.4
390   65  CC   4.594647   H116A  2023-03-03    H1  16    A    4  10       8.0
391   65  CC   4.594647   H234B  2023-03-03    H2  34    B    4   7       3.5
392   65  CC   4.594647   H236C  2023-03-03    H2  36    C    2   1       5.9
393   65  CC   4.594647   J112B  2023-03-03    J1  12    B    4   4       5.8
394   65  CC   4.594647   J228A  2023-03-03    J2  28    A    5   6       4.3
395   65  CC   4.594647   J236C  2023-03-03    J2  36    C    5   5       5.2
396   65  CC   4.594647   J239C  2023-03-03    J2  39    C    2   3       8.2
397   65  CC   4.594647    L18A  2023-03-03    L1   8    A    4   6       5.5
398   65  CC   4.594647   L232A  2023-03-03    L2  32    A    5   3       6.3
399   65  CC   4.594647   L117C  2023-03-03    L1  17    C    7   3       4.8
400   65  BH   4.400758    D18B  2023-03-03    D1   8    B    6   7       3.3
401   65  BH   4.400758   D120A  2023-03-03    D1  20    A    4   2       6.5
402   65  BH   4.400758   D226B  2023-03-03    D2  26    B    5  15       4.2
403   65  BH   4.400758   D230B  2023-03-03    D2  30    B    4   3       4.9
404   65  BH   4.400758   D233B  2023-03-03    D2  33    B    6   8       4.0
405   65  BH   4.400758    D16C  2023-03-03    D1   6    C    4  12       5.9
406   65  BH   4.400758    D19C  2023-03-03    D1   9    C    1   8       2.7
407   65  BH   4.400758   D119D  2023-03-03    D1  19    D    5   9       1.4
408   65  BH   4.400758   D229C  2023-03-03    D2  29    C    4  13       4.7
409   65  BH   4.400758   D229D  2023-03-03    D2  29    D    2   9       4.8
410   65  BH   4.400758   D235D  2023-03-03    D2  35    D    3   6       5.2
411   65  BH   4.400758   F115A  2023-03-03    F1  15    A    6   9       5.5
412   65  BH   4.400758   F229B  2023-03-03    F2  29    B    4   5       5.0
413   65  BH   4.400758   F240B  2023-03-03    F2  40    B    5  12       4.3
414   65  BH   4.400758    F15C  2023-03-03    F1   5    C    5  26       2.9
415   65  BH   4.400758   F113C  2023-03-03    F1  13    C    3   7       3.8
416   65  BH   4.400758   F115C  2023-03-03    F1  15    C    4   4       4.9
417   65  BH   4.400758   F235C  2023-03-03    F2  35    C    5  27       4.2
418   65  BH   4.400758   F235D  2023-03-03    F2  35    D    2   1       4.5
419   65  BH   4.400758    H14B  2023-03-03    H1   4    B    2  12       5.1
420   65  BH   4.400758    H18B  2023-03-03    H1   8    B    3   9       3.5
421   65  BH   4.400758   H114A  2023-03-03    H1  14    A    5   3       6.9
422   65  BH   4.400758   H117A  2023-03-03    H1  17    A    4   6       3.3
423   65  BH   4.400758   H118B  2023-03-03    H1  18    B    5  18       5.0
424   65  BH   4.400758   H224B  2023-03-03    H2  24    B    5  29       9.3
425   65  BH   4.400758   H227A  2023-03-03    H2  27    A    5   4       3.9
426   65  BH   4.400758   H229B  2023-03-03    H2  29    B    5  19       4.4
427   65  BH   4.400758   H233B  2023-03-03    H2  33    B    4   7       2.5
428   65  BH   4.400758   H239B  2023-03-03    H2  39    B    6   2       3.5
429   65  BH   4.400758    H15C  2023-03-03    H1   5    C    5  28       5.3
430   65  BH   4.400758   H110C  2023-03-03    H1  10    C    6  11       3.7
431   65  BH   4.400758   H115D  2023-03-03    H1  15    D    5  13       4.8
432   65  BH   4.400758   H119C  2023-03-03    H1  19    C    7  21       4.6
433   65  BH   4.400758   H222D  2023-03-03    H2  22    D    2   3       3.1
434   65  BH   4.400758   H223C  2023-03-03    H2  23    C    7  22       3.4
435   65  BH   4.400758   H235C  2023-03-03    H2  35    C    2  13       2.8
436   65  BH   4.400758    J14B  2023-03-03    J1   4    B    2   4       3.6
437   65  BH   4.400758   J119B  2023-03-03    J1  19    B    3  11       5.6
438   65  BH   4.400758   J233A  2023-03-03    J2  33    A    3   2       4.6
439   65  BH   4.400758   J237B  2023-03-03    J2  37    B    1   7       4.8
440   65  BH   4.400758    J19C  2023-03-03    J1   9    C    5  30       5.6
441   65  BH   4.400758   J110C  2023-03-03    J1  10    C    3   1       2.2
442   65  BH   4.400758   J114D  2023-03-03    J1  14    D    4   8       3.5
443   65  BH   4.400758   J119C  2023-03-03    J1  19    C    5  20       5.2
444   65  BH   4.400758   J223D  2023-03-03    J2  23    D    4   9       4.9
445   65  BH   4.400758   J225D  2023-03-03    J2  25    D    2   5       1.9
446   65  BH   4.400758   J230C  2023-03-03    J2  30    C    4  14       6.3
447   65  BH   4.400758   J231C  2023-03-03    J2  31    C    5   6       7.8
448   65  BH   4.400758   J231D  2023-03-03    J2  31    D    7  14       7.7
449   65  BH   4.400758   J232D  2023-03-03    J2  32    D    3  12       3.5
450   65  BH   4.400758    L15B  2023-03-03    L1   5    B    2   6       3.8
451   65  BH   4.400758    L18B  2023-03-03    L1   8    B    7   2       2.3
452   65  BH   4.400758   L110A  2023-03-03    L1  10    A    2  14       1.2
453   65  BH   4.400758   L111B  2023-03-03    L1  11    B    3   3       2.5
454   65  BH   4.400758    L14D  2023-03-03    L1   4    D    4  10       3.6
455   65  BH   4.400758   L113C  2023-03-03    L1  13    C    3  13       2.4
456   65  BH   4.400758   L113D  2023-03-03    L1  13    D    5  22       4.6
457   65  BH   4.400758   L116C  2023-03-03    L1  16    C    6   5       4.8
458   65  BH   4.400758   L228C  2023-03-03    L2  28    C    6   6       2.2
459   65  BH   4.400758   L228D  2023-03-03    L2  28    D    2   7       3.5
460   65  BH   4.400758   L229D  2023-03-03    L2  29    D    3  14       2.1
461   65  BH   4.400758   L230C  2023-03-03    L2  30    C    4  11       5.1
462   65  BH   4.400758   L234C  2023-03-03    L2  34    C    5   8       4.4
463   65  BH   4.400758   L240D  2023-03-03    L2  40    D    7  13       4.3
464   65  BH   4.400758   L241D  2023-03-03    L2  41    D    4   1       3.2
465   65 DPR   4.903883   D225A  2023-03-03    D2  25    A    2   7       9.3
466   65 DPR   4.903883   D227B  2023-03-03    D2  27    B    4   9       7.9
467   65 DPR   4.903883   F232A  2023-03-03    F2  32    A    6   2       2.7
468   65 DPR   4.903883    F16C  2023-03-03    F1   6    C    6   9       4.3
469   65 DPR   4.903883   F240C  2023-03-03    F2  40    C    7   4       5.8
470   65 DPR   4.903883   H117B  2023-03-03    H1  17    B    5   2      10.1
471   65 DPR   4.903883    H17C  2023-03-03    H1   7    C    2  10       9.2
472   65 DPR   4.903883   H225C  2023-03-03    H2  25    C    2   6       7.6
473   65 DPR   4.903883    J16B  2023-03-03    J1   6    B    6  10       4.5
474   65 DPR   4.903883   J232B  2023-03-03    J2  32    B    6  11       4.4
475   65 DPR   4.903883   L237B  2023-03-03    L2  37    B    6   3       5.1
476   66  CC   4.349152    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
477   66  BH   4.348820    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
478   66 DPR   4.846899    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
479   67  CC   4.468296    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
480   67  BH   4.639043    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
481   67 DPR   4.780840    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
482   68  CC   4.639013    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
483   68  BH   4.680593    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
484   68 DPR   5.183934    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
485   69  CC   4.612849    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
486   69  BH   4.496744    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
487   69 DPR   4.792970    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
488   70  CC   4.898392    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
489   70  BH   4.505360    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
490   70 DPR   5.156766    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
491   71  CC   4.937836    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
492   71  BH   4.752797    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
493   71 DPR   5.070183    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
494   72  CC   4.546545    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
495   72  BH   4.977534    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
496   72 DPR   5.201028    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
497   73  CC   4.909303    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
498   73  BH   4.365176    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
499   73 DPR   5.325978    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
500   74  CC   4.884729    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
501   74  BH   4.884963    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
502   74 DPR   5.325565    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
503   75  CC   5.185708    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
504   75  BH   4.746108    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
505   75 DPR   5.717116    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
506   76  CC   5.282582    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
507   76  BH   5.057293    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
508   76 DPR   5.482464    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
509   77  CC   5.234717    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
510   77  BH   5.285357    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
511   77 DPR   5.892511    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
512   78  CC   5.200374    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
513   78  BH   5.329238    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
514   78 DPR   5.957380    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
515   79  CC   5.451101   D112A  2023-03-17    D1  12    A    7   4       5.4
516   79  CC   5.451101   D239A  2023-03-17    D2  39    A    5   7       3.9
517   79  CC   5.451101   D239B  2023-03-17    D2  39    B    4   5       5.2
518   79  CC   5.451101   F117D  2023-03-17    F1  17    D    3  13       7.2
519   79  CC   5.451101   H111B  2023-03-17    H1  11    B    7   2      12.3
520   79  CC   5.451101   H115B  2023-03-17    H1  15    B    3   3       9.7
521   79  CC   5.451101   H116A  2023-03-17    H1  16    A    4  10       9.1
522   79  CC   5.451101   H234B  2023-03-17    H2  34    B    4   7       3.4
523   79  CC   5.451101   H236C  2023-03-17    H2  36    C    2   1       6.1
524   79  CC   5.451101   J112B  2023-03-17    J1  12    B    4   4       6.3
525   79  CC   5.451101   J236C  2023-03-17    J2  36    C    5   5       7.2
526   79  CC   5.451101    L18A  2023-03-17    L1   8    A    4   6       4.2
527   79  CC   5.451101   L232A  2023-03-17    L2  32    A    5   3       5.9
528   79  CC   5.451101   L117C  2023-03-17    L1  17    C    7   3       5.4
529   79  BH   5.180753   D226B  2023-03-17    D2  26    B    5  15       5.3
530   79  BH   5.180753   D230B  2023-03-17    D2  30    B    4   3       6.1
531   79  BH   5.180753   D233B  2023-03-17    D2  33    B    6   8       5.6
532   79  BH   5.180753    D19C  2023-03-17    D1   9    C    1   8       2.6
533   79  BH   5.180753   D119D  2023-03-17    D1  19    D    5   9       2.1
534   79  BH   5.180753   D229C  2023-03-17    D2  29    C    4  13       4.3
535   79  BH   5.180753   D229D  2023-03-17    D2  29    D    2   9       4.9
536   79  BH   5.180753   D235D  2023-03-17    D2  35    D    3   6       5.3
537   79  BH   5.180753   F115A  2023-03-17    F1  15    A    6   9       8.0
538   79  BH   5.180753   F229B  2023-03-17    F2  29    B    4   5       5.9
539   79  BH   5.180753   F240B  2023-03-17    F2  40    B    5  12       5.1
540   79  BH   5.180753    F15C  2023-03-17    F1   5    C    5  26       3.0
541   79  BH   5.180753   F113C  2023-03-17    F1  13    C    3   7       4.1
542   79  BH   5.180753   F115C  2023-03-17    F1  15    C    4   4       4.5
543   79  BH   5.180753   F235C  2023-03-17    F2  35    C    5  27       5.5
544   79  BH   5.180753   F235D  2023-03-17    F2  35    D    2   1       4.6
545   79  BH   5.180753    H14B  2023-03-17    H1   4    B    2  12       5.5
546   79  BH   5.180753    H18B  2023-03-17    H1   8    B    3   9       6.5
547   79  BH   5.180753   H114A  2023-03-17    H1  14    A    5   3       6.2
548   79  BH   5.180753   H117A  2023-03-17    H1  17    A    4   6       3.6
549   79  BH   5.180753   H118B  2023-03-17    H1  18    B    5  18       5.0
550   79  BH   5.180753   H224B  2023-03-17    H2  24    B    5  29       8.5
551   79  BH   5.180753   H233B  2023-03-17    H2  33    B    4   7       3.5
552   79  BH   5.180753   H239B  2023-03-17    H2  39    B    6   2       4.0
553   79  BH   5.180753   H110C  2023-03-17    H1  10    C    6  11       5.6
554   79  BH   5.180753   H115D  2023-03-17    H1  15    D    5  13       5.6
555   79  BH   5.180753   H119C  2023-03-17    H1  19    C    7  21       6.5
556   79  BH   5.180753   H222D  2023-03-17    H2  22    D    2   3       4.4
557   79  BH   5.180753   H223C  2023-03-17    H2  23    C    7  22       4.9
558   79  BH   5.180753   H235C  2023-03-17    H2  35    C    2  13       4.8
559   79  BH   5.180753    J14B  2023-03-17    J1   4    B    2   4       4.5
560   79  BH   5.180753   J119B  2023-03-17    J1  19    B    3  11       5.3
561   79  BH   5.180753   J233A  2023-03-17    J2  33    A    3   2       4.7
562   79  BH   5.180753   J237B  2023-03-17    J2  37    B    1   7       5.3
563   79  BH   5.180753    J19C  2023-03-17    J1   9    C    5  30       8.9
564   79  BH   5.180753   J110C  2023-03-17    J1  10    C    3   1       3.9
565   79  BH   5.180753   J114D  2023-03-17    J1  14    D    4   8       3.8
566   79  BH   5.180753   J119C  2023-03-17    J1  19    C    5  20       7.1
567   79  BH   5.180753   J223D  2023-03-17    J2  23    D    4   9       5.2
568   79  BH   5.180753   J225D  2023-03-17    J2  25    D    2   5       3.0
569   79  BH   5.180753   J230C  2023-03-17    J2  30    C    4  14       8.3
570   79  BH   5.180753   J231C  2023-03-17    J2  31    C    5   6       9.2
571   79  BH   5.180753   J232D  2023-03-17    J2  32    D    3  12       4.5
572   79  BH   5.180753    L15B  2023-03-17    L1   5    B    2   6       4.2
573   79  BH   5.180753   L110A  2023-03-17    L1  10    A    2  14       2.0
574   79  BH   5.180753    L14D  2023-03-17    L1   4    D    4  10       3.8
575   79  BH   5.180753   L113C  2023-03-17    L1  13    C    3  13       3.6
576   79  BH   5.180753   L113D  2023-03-17    L1  13    D    5  22       6.5
577   79  BH   5.180753   L116C  2023-03-17    L1  16    C    6   5       7.9
578   79  BH   5.180753   L228D  2023-03-17    L2  28    D    2   7       5.1
579   79  BH   5.180753   L229D  2023-03-17    L2  29    D    3  14       3.1
580   79  BH   5.180753   L230C  2023-03-17    L2  30    C    4  11       6.6
581   79  BH   5.180753   L234C  2023-03-17    L2  34    C    5   8       4.9
582   79  BH   5.180753   L240D  2023-03-17    L2  40    D    7  13       6.0
583   79  BH   5.180753   L241D  2023-03-17    L2  41    D    4   1       4.6
584   79 DPR   6.359643   F232A  2023-03-17    F2  32    A    6   2       2.5
585   79 DPR   6.359643    F16C  2023-03-17    F1   6    C    6   9       6.6
586   79 DPR   6.359643   F240C  2023-03-17    F2  40    C    7   4      10.2
587   79 DPR   6.359643   H117B  2023-03-17    H1  17    B    5   2      17.0
588   79 DPR   6.359643    H17C  2023-03-17    H1   7    C    2  10      11.4
589   79 DPR   6.359643   H225C  2023-03-17    H2  25    C    2   6      12.4
590   79 DPR   6.359643    J16B  2023-03-17    J1   6    B    6  10       7.4
591   79 DPR   6.359643   J232B  2023-03-17    J2  32    B    6  11       4.2
592   79 DPR   6.359643   L237B  2023-03-17    L2  37    B    6   3       6.6
593   80  CC   5.458317    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
594   80  BH   5.346463    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
595   80 DPR   6.571951    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
596   81  CC   5.505370    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
597   81  BH   5.132736    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
598   81 DPR   6.589102    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
599   82  CC   6.032080    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
600   82  BH   5.471146    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
601   82 DPR   6.846597    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
602   83  CC   6.179954    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
603   83  BH   5.445137    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
604   83 DPR   6.849081    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
605   84  CC   6.098192    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
606   84  BH   5.677515    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
607   84 DPR   6.840962    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
608   85  CC   6.364802    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
609   85  BH   6.009689    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
610   85 DPR   7.424297    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
611   86  CC   6.497986   D112A  2023-03-24    D1  12    A    7   4       5.5
612   86  CC   6.497986   H111B  2023-03-24    H1  11    B    7   2      12.6
613   86  CC   6.497986   H115B  2023-03-24    H1  15    B    3   3      12.5
614   86  CC   6.497986   H116A  2023-03-24    H1  16    A    4  10      10.1
615   86  CC   6.497986   H234B  2023-03-24    H2  34    B    4   7       3.1
616   86  CC   6.497986   H236C  2023-03-24    H2  36    C    2   1       6.0
617   86  CC   6.497986   J112B  2023-03-24    J1  12    B    4   4       9.8
618   86  CC   6.497986   J236C  2023-03-24    J2  36    C    5   5       7.0
619   86  CC   6.497986    L18A  2023-03-24    L1   8    A    4   6       4.6
620   86  CC   6.497986   L232A  2023-03-24    L2  32    A    5   3       6.3
621   86  CC   6.497986   L117C  2023-03-24    L1  17    C    7   3       5.5
622   86  BH   6.327779   D226B  2023-03-24    D2  26    B    5  15       5.4
623   86  BH   6.327779   D230B  2023-03-24    D2  30    B    4   3       3.3
624   86  BH   6.327779   D233B  2023-03-24    D2  33    B    6   8       3.9
625   86  BH   6.327779   D119D  2023-03-24    D1  19    D    5   9       3.5
626   86  BH   6.327779   D229C  2023-03-24    D2  29    C    4  13       4.8
627   86  BH   6.327779   D229D  2023-03-24    D2  29    D    2   9       6.1
628   86  BH   6.327779   D235D  2023-03-24    D2  35    D    3   6       5.7
629   86  BH   6.327779   F115A  2023-03-24    F1  15    A    6   9       4.5
630   86  BH   6.327779   F229B  2023-03-24    F2  29    B    4   5       5.6
631   86  BH   6.327779   F240B  2023-03-24    F2  40    B    5  12       5.7
632   86  BH   6.327779   F115C  2023-03-24    F1  15    C    4   4       4.9
633   86  BH   6.327779   F235D  2023-03-24    F2  35    D    2   1       5.1
634   86  BH   6.327779    H14B  2023-03-24    H1   4    B    2  12       7.3
635   86  BH   6.327779    H18B  2023-03-24    H1   8    B    3   9       5.4
636   86  BH   6.327779   H114A  2023-03-24    H1  14    A    5   3       7.1
637   86  BH   6.327779   H224B  2023-03-24    H2  24    B    5  29       9.8
638   86  BH   6.327779   H233B  2023-03-24    H2  33    B    4   7       4.5
639   86  BH   6.327779   H110C  2023-03-24    H1  10    C    6  11       8.0
640   86  BH   6.327779   H115D  2023-03-24    H1  15    D    5  13       4.9
641   86  BH   6.327779   H119C  2023-03-24    H1  19    C    7  21       7.0
642   86  BH   6.327779   H222D  2023-03-24    H2  22    D    2   3       5.5
643   86  BH   6.327779   H223C  2023-03-24    H2  23    C    7  22       7.4
644   86  BH   6.327779   H235C  2023-03-24    H2  35    C    2  13       3.8
645   86  BH   6.327779    J14B  2023-03-24    J1   4    B    2   4       3.9
646   86  BH   6.327779   J119B  2023-03-24    J1  19    B    3  11       6.4
647   86  BH   6.327779   J233A  2023-03-24    J2  33    A    3   2       5.9
648   86  BH   6.327779   J237B  2023-03-24    J2  37    B    1   7       7.2
649   86  BH   6.327779    J19C  2023-03-24    J1   9    C    5  30       9.7
650   86  BH   6.327779   J110C  2023-03-24    J1  10    C    3   1       4.5
651   86  BH   6.327779   J119C  2023-03-24    J1  19    C    5  20       7.5
652   86  BH   6.327779   J223D  2023-03-24    J2  23    D    4   9       5.5
653   86  BH   6.327779   J225D  2023-03-24    J2  25    D    2   5       3.0
654   86  BH   6.327779   J230C  2023-03-24    J2  30    C    4  14       8.5
655   86  BH   6.327779   J231D  2023-03-24    J2  31    D    7  14       8.3
656   86  BH   6.327779   J232D  2023-03-24    J2  32    D    3  12       4.8
657   86  BH   6.327779    L15B  2023-03-24    L1   5    B    2   6       4.5
658   86  BH   6.327779   L110A  2023-03-24    L1  10    A    2  14       3.1
659   86  BH   6.327779    L14D  2023-03-24    L1   4    D    4  10       4.1
660   86  BH   6.327779   L113C  2023-03-24    L1  13    C    3  13       2.8
661   86  BH   6.327779   L113D  2023-03-24    L1  13    D    5  22       6.5
662   86  BH   6.327779   L228D  2023-03-24    L2  28    D    2   7       5.5
663   86  BH   6.327779   L229D  2023-03-24    L2  29    D    3  14       3.3
664   86  BH   6.327779   L230C  2023-03-24    L2  30    C    4  11       5.1
665   86  BH   6.327779   L240D  2023-03-24    L2  40    D    7  13       5.3
666   86  BH   6.327779   L241D  2023-03-24    L2  41    D    4   1       4.7
667   86 DPR   7.657962    F16C  2023-03-24    F1   6    C    6   9       9.1
668   86 DPR   7.657962   F240C  2023-03-24    F2  40    C    7   4      15.4
669   86 DPR   7.657962   H117B  2023-03-24    H1  17    B    5   2      21.7
670   86 DPR   7.657962    H17C  2023-03-24    H1   7    C    2  10      12.5
671   86 DPR   7.657962   H225C  2023-03-24    H2  25    C    2   6      18.4
672   86 DPR   7.657962    J16B  2023-03-24    J1   6    B    6  10      10.8
673   86 DPR   7.657962   J232B  2023-03-24    J2  32    B    6  11       4.4
674   86 DPR   7.657962   L237B  2023-03-24    L2  37    B    6   3       6.4
675   87  CC   6.800994    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
676   87  BH   6.141755    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
677   87 DPR   7.942002    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
678   88  CC   7.056536    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
679   88  BH   6.390540    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
680   88 DPR   8.180654    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
681   89  CC   7.436471    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
682   89  BH   6.416927    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
683   89 DPR   8.639910    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
684   90  CC   7.226212    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
685   90  BH   6.956599    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
686   90 DPR   8.835446    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
687   91  CC   7.191260    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
688   91  BH   6.967425    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
689   91 DPR   9.470547    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
690   92  CC   7.642828    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
691   92  BH   6.829581    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
692   92 DPR   9.649409    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
693   93  CC   7.902681    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
694   93  BH   7.174991    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
695   93 DPR   9.897973    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
696   94  CC   8.341474    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
697   94  BH   7.405078    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
698   94 DPR  10.330867    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
699   95  CC   8.262334    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
700   95  BH   7.509571    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
701   95 DPR  10.427684    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
702   96  CC   8.740141   H111B  2023-04-03    H1  11    B    7   2      12.1
703   96  CC   8.740141   H115B  2023-04-03    H1  15    B    3   3      14.9
704   96  CC   8.740141   H116A  2023-04-03    H1  16    A    4  10      10.2
705   96  CC   8.740141   J112B  2023-04-03    J1  12    B    4   4       7.5
706   96  CC   8.740141   J236C  2023-04-03    J2  36    C    5   5       8.3
707   96  CC   8.740141   L117C  2023-04-03    L1  17    C    7   3       6.2
708   96  BH   7.977976   D226B  2023-04-03    D2  26    B    5  15       4.9
709   96  BH   7.977976   D230B  2023-04-03    D2  30    B    4   3       3.8
710   96  BH   7.977976   D119D  2023-04-03    D1  19    D    5   9       3.2
711   96  BH   7.977976   D229C  2023-04-03    D2  29    C    4  13       4.9
712   96  BH   7.977976   D229D  2023-04-03    D2  29    D    2   9       4.1
713   96  BH   7.977976   D235D  2023-04-03    D2  35    D    3   6       6.6
714   96  BH   7.977976   F115A  2023-04-03    F1  15    A    6   9       9.9
715   96  BH   7.977976   F229B  2023-04-03    F2  29    B    4   5       5.9
716   96  BH   7.977976   F240B  2023-04-03    F2  40    B    5  12       5.6
717   96  BH   7.977976   F115C  2023-04-03    F1  15    C    4   4       5.9
718   96  BH   7.977976   F235D  2023-04-03    F2  35    D    2   1       5.5
719   96  BH   7.977976    H14B  2023-04-03    H1   4    B    2  12       0.7
720   96  BH   7.977976    H18B  2023-04-03    H1   8    B    3   9       7.1
721   96  BH   7.977976   H114A  2023-04-03    H1  14    A    5   3       8.4
722   96  BH   7.977976   H224B  2023-04-03    H2  24    B    5  29       9.4
723   96  BH   7.977976   H233B  2023-04-03    H2  33    B    4   7       4.8
724   96  BH   7.977976   H110C  2023-04-03    H1  10    C    6  11      10.6
725   96  BH   7.977976   H119C  2023-04-03    H1  19    C    7  21       6.5
726   96  BH   7.977976   H223C  2023-04-03    H2  23    C    7  22       8.4
727   96  BH   7.977976   H235C  2023-04-03    H2  35    C    2  13       3.8
728   96  BH   7.977976    J14B  2023-04-03    J1   4    B    2   4       5.1
729   96  BH   7.977976   J119B  2023-04-03    J1  19    B    3  11       5.9
730   96  BH   7.977976   J233A  2023-04-03    J2  33    A    3   2       5.8
731   96  BH   7.977976   J237B  2023-04-03    J2  37    B    1   7       4.4
732   96  BH   7.977976    J19C  2023-04-03    J1   9    C    5  30      10.2
733   96  BH   7.977976   J110C  2023-04-03    J1  10    C    3   1       5.1
734   96  BH   7.977976   J119C  2023-04-03    J1  19    C    5  20      10.1
735   96  BH   7.977976   J223D  2023-04-03    J2  23    D    4   9       5.7
736   96  BH   7.977976   J225D  2023-04-03    J2  25    D    2   5       6.3
737   96  BH   7.977976   J230C  2023-04-03    J2  30    C    4  14       9.4
738   96  BH   7.977976   J231D  2023-04-03    J2  31    D    7  14      12.1
739   96  BH   7.977976   J232D  2023-04-03    J2  32    D    3  12       5.2
740   96  BH   7.977976    L15B  2023-04-03    L1   5    B    2   6       3.9
741   96  BH   7.977976   L110A  2023-04-03    L1  10    A    2  14       3.6
742   96  BH   7.977976    L14D  2023-04-03    L1   4    D    4  10       7.1
743   96  BH   7.977976   L113C  2023-04-03    L1  13    C    3  13       3.6
744   96  BH   7.977976   L113D  2023-04-03    L1  13    D    5  22       7.3
745   96  BH   7.977976   L228D  2023-04-03    L2  28    D    2   7       7.1
746   96  BH   7.977976   L229D  2023-04-03    L2  29    D    3  14       4.1
747   96  BH   7.977976   L230C  2023-04-03    L2  30    C    4  11       4.3
748   96  BH   7.977976   L240D  2023-04-03    L2  40    D    7  13       6.6
749   96  BH   7.977976   L241D  2023-04-03    L2  41    D    4   1       5.8
750   96 DPR  11.437417    F16C  2023-04-03    F1   6    C    6   9      15.4
751   96 DPR  11.437417   F240C  2023-04-03    F2  40    C    7   4      22.5
752   96 DPR  11.437417   H117B  2023-04-03    H1  17    B    5   2      25.9
753   96 DPR  11.437417    H17C  2023-04-03    H1   7    C    2  10      13.3
754   96 DPR  11.437417   H225C  2023-04-03    H2  25    C    2   6      27.7
755   97  CC   8.915929    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
756   97  BH   8.189641    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
757   97 DPR  11.612699    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
758   98  CC   9.221648    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
759   98  BH   8.388396    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
760   98 DPR  12.030501    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
761   99  CC   9.678219    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
762   99  BH   8.668453    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
763   99 DPR  12.364770    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
764  100  CC  10.073818    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
765  100  BH   8.943405    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
766  100 DPR  12.569543    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
767  101  CC  10.167191    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
768  101  BH   9.408457    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
769  101 DPR  13.699703    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
770  102  CC  10.574377    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
771  102  BH   9.381471    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
772  102 DPR  13.720811    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
773  103  CC  10.835829   H111B  2023-04-10    H1  11    B    7   2      16.2
774  103  CC  10.835829   H115B  2023-04-10    H1  15    B    3   3      23.8
775  103  CC  10.835829   H116A  2023-04-10    H1  16    A    4  10      12.3
776  103  CC  10.835829   J112B  2023-04-10    J1  12    B    4   4      10.1
777  103  CC  10.835829   J236C  2023-04-10    J2  36    C    5   5       5.9
778  103  BH   9.927091   D226B  2023-04-10    D2  26    B    5  15       5.7
779  103  BH   9.927091   D230B  2023-04-10    D2  30    B    4   3       6.8
780  103  BH   9.927091   D119D  2023-04-10    D1  19    D    5   9       3.6
781  103  BH   9.927091   D229C  2023-04-10    D2  29    C    4  13       5.8
782  103  BH   9.927091   D229D  2023-04-10    D2  29    D    2   9       4.3
783  103  BH   9.927091   D235D  2023-04-10    D2  35    D    3   6       8.4
784  103  BH   9.927091   F115A  2023-04-10    F1  15    A    6   9      17.4
785  103  BH   9.927091   F229B  2023-04-10    F2  29    B    4   5       7.2
786  103  BH   9.927091   F240B  2023-04-10    F2  40    B    5  12       6.6
787  103  BH   9.927091   F115C  2023-04-10    F1  15    C    4   4       5.5
788  103  BH   9.927091   F235D  2023-04-10    F2  35    D    2   1       5.7
789  103  BH   9.927091    H14B  2023-04-10    H1   4    B    2  12       6.1
790  103  BH   9.927091    H18B  2023-04-10    H1   8    B    3   9       8.8
791  103  BH   9.927091   H114A  2023-04-10    H1  14    A    5   3      12.3
792  103  BH   9.927091   H224B  2023-04-10    H2  24    B    5  29      13.4
793  103  BH   9.927091   H233B  2023-04-10    H2  33    B    4   7       5.7
794  103  BH   9.927091   H110C  2023-04-10    H1  10    C    6  11      13.2
795  103  BH   9.927091   H119C  2023-04-10    H1  19    C    7  21       5.1
796  103  BH   9.927091   H223C  2023-04-10    H2  23    C    7  22       8.6
797  103  BH   9.927091   H235C  2023-04-10    H2  35    C    2  13       3.9
798  103  BH   9.927091    J14B  2023-04-10    J1   4    B    2   4       6.0
799  103  BH   9.927091   J119B  2023-04-10    J1  19    B    3  11       6.9
800  103  BH   9.927091   J233A  2023-04-10    J2  33    A    3   2       6.4
801  103  BH   9.927091   J237B  2023-04-10    J2  37    B    1   7      12.2
802  103  BH   9.927091    J19C  2023-04-10    J1   9    C    5  30       8.5
803  103  BH   9.927091   J110C  2023-04-10    J1  10    C    3   1       7.0
804  103  BH   9.927091   J119C  2023-04-10    J1  19    C    5  20      11.4
805  103  BH   9.927091   J225D  2023-04-10    J2  25    D    2   5       7.7
806  103  BH   9.927091   J230C  2023-04-10    J2  30    C    4  14       9.0
807  103  BH   9.927091   J231D  2023-04-10    J2  31    D    7  14      10.2
808  103  BH   9.927091   J232D  2023-04-10    J2  32    D    3  12       6.0
809  103  BH   9.927091   L110A  2023-04-10    L1  10    A    2  14       2.1
810  103  BH   9.927091    L14D  2023-04-10    L1   4    D    4  10       4.1
811  103  BH   9.927091   L113C  2023-04-10    L1  13    C    3  13       2.5
812  103  BH   9.927091   L113D  2023-04-10    L1  13    D    5  22       9.0
813  103  BH   9.927091   L228D  2023-04-10    L2  28    D    2   7       6.3
814  103  BH   9.927091   L229D  2023-04-10    L2  29    D    3  14       3.6
815  103  BH   9.927091   L240D  2023-04-10    L2  40    D    7  13       8.4
816  103  BH   9.927091   L241D  2023-04-10    L2  41    D    4   1       4.4
817  103 DPR  14.824650    F16C  2023-04-10    F1   6    C    6   9      16.6
818  103 DPR  14.824650   H117B  2023-04-10    H1  17    B    5   2      38.2
819  103 DPR  14.824650   H225C  2023-04-10    H2  25    C    2   6      37.4
820  103 DPR  14.824650    J16B  2023-04-10    J1   6    B    6  10      23.6
821  104  CC  11.204196    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
822  104  BH  10.306400    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
823  104 DPR  15.058698    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
824  105  CC  12.088986    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
825  105  BH  10.787903    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
826  105 DPR  16.011845    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
827  106  CC  12.060390    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
828  106  BH  10.545755    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
829  106 DPR  16.603843    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
830  107  CC  12.782610    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
831  107  BH  11.437026    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
832  107 DPR  16.950690    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
833  108  CC  13.382431    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
834  108  BH  11.748042    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
835  108 DPR  17.728514    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
836  109  CC  13.550854    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
837  109  BH  12.023989    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
838  109 DPR  18.607704    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
839  110  CC  14.048636   H111B  2023-04-17    H1  11    B    7   2      14.4
840  110  CC  14.048636   H115B  2023-04-17    H1  15    B    3   3      32.6
841  110  CC  14.048636   H116A  2023-04-17    H1  16    A    4  10      15.1
842  110  CC  14.048636   J112B  2023-04-17    J1  12    B    4   4       9.0
843  110  BH  12.724186   D226B  2023-04-17    D2  26    B    5  15       7.9
844  110  BH  12.724186   D230B  2023-04-17    D2  30    B    4   3       8.1
845  110  BH  12.724186   D119D  2023-04-17    D1  19    D    5   9       3.3
846  110  BH  12.724186   D229C  2023-04-17    D2  29    C    4  13       6.8
847  110  BH  12.724186   D229D  2023-04-17    D2  29    D    2   9       7.1
848  110  BH  12.724186   D235D  2023-04-17    D2  35    D    3   6      10.2
849  110  BH  12.724186   F115A  2023-04-17    F1  15    A    6   9      22.9
850  110  BH  12.724186   F229B  2023-04-17    F2  29    B    4   5       8.2
851  110  BH  12.724186   F240B  2023-04-17    F2  40    B    5  12       8.0
852  110  BH  12.724186   F235D  2023-04-17    F2  35    D    2   1       7.3
853  110  BH  12.724186    H14B  2023-04-17    H1   4    B    2  12      11.4
854  110  BH  12.724186    H18B  2023-04-17    H1   8    B    3   9      15.2
855  110  BH  12.724186   H114A  2023-04-17    H1  14    A    5   3      13.1
856  110  BH  12.724186   H224B  2023-04-17    H2  24    B    5  29      19.2
857  110  BH  12.724186   H233B  2023-04-17    H2  33    B    4   7       6.2
858  110  BH  12.724186   H110C  2023-04-17    H1  10    C    6  11      18.3
859  110  BH  12.724186   H223C  2023-04-17    H2  23    C    7  22      11.8
860  110  BH  12.724186    J14B  2023-04-17    J1   4    B    2   4      12.9
861  110  BH  12.724186   J119B  2023-04-17    J1  19    B    3  11      12.0
862  110  BH  12.724186   J233A  2023-04-17    J2  33    A    3   2       9.0
863  110  BH  12.724186   J237B  2023-04-17    J2  37    B    1   7      18.4
864  110  BH  12.724186    J19C  2023-04-17    J1   9    C    5  30       9.8
865  110  BH  12.724186   J110C  2023-04-17    J1  10    C    3   1      13.4
866  110  BH  12.724186   J119C  2023-04-17    J1  19    C    5  20      18.4
867  110  BH  12.724186   J223D  2023-04-17    J2  23    D    4   9       6.1
868  110  BH  12.724186   J225D  2023-04-17    J2  25    D    2   5       8.8
869  110  BH  12.724186   J230C  2023-04-17    J2  30    C    4  14      13.1
870  110  BH  12.724186   J231D  2023-04-17    J2  31    D    7  14      11.2
871  110  BH  12.724186   J232D  2023-04-17    J2  32    D    3  12      10.8
872  110  BH  12.724186   L110A  2023-04-17    L1  10    A    2  14      11.1
873  110  BH  12.724186    L14D  2023-04-17    L1   4    D    4  10      14.4
874  110  BH  12.724186   L113C  2023-04-17    L1  13    C    3  13      16.0
875  110  BH  12.724186   L113D  2023-04-17    L1  13    D    5  22      12.4
876  110  BH  12.724186   L228D  2023-04-17    L2  28    D    2   7       7.1
877  110  BH  12.724186   L229D  2023-04-17    L2  29    D    3  14       5.6
878  110  BH  12.724186   L240D  2023-04-17    L2  40    D    7  13      11.8
879  110  BH  12.724186   L241D  2023-04-17    L2  41    D    4   1       6.9
880  110 DPR  19.633113    F16C  2023-04-17    F1   6    C    6   9      16.1
881  110 DPR  19.633113   H117B  2023-04-17    H1  17    B    5   2      42.7
882  110 DPR  19.633113   H225C  2023-04-17    H2  25    C    2   6      36.2
883  110 DPR  19.633113    J16B  2023-04-17    J1   6    B    6  10      24.5
884  111  CC  14.261759    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
885  111  BH  13.206855    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
886  111 DPR  20.006144    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
887  112  CC  14.639875    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
888  112  BH  13.614070    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
889  112 DPR  20.920895    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
890  113  CC  15.838977    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
891  113  BH  13.809135    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
892  113 DPR  21.976439    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
893  114  CC  16.328254    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
894  114  BH  14.360145    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
895  114 DPR  22.987894    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
896  115  CC  16.567823    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
897  115  BH  14.935287    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
898  115 DPR  23.519926    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
899  116  CC  17.799015    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
900  116  BH  15.542638    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
901  116 DPR  24.767556    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
902  117  CC  17.652498   H111B  2023-04-24    H1  11    B    7   2      27.1
903  117  CC  17.652498   H115B  2023-04-24    H1  15    B    3   3      37.0
904  117  CC  17.652498   H116A  2023-04-24    H1  16    A    4  10      15.9
905  117  BH  15.854418   D226B  2023-04-24    D2  26    B    5  15      15.9
906  117  BH  15.854418   D230B  2023-04-24    D2  30    B    4   3      12.6
907  117  BH  15.854418   D119D  2023-04-24    D1  19    D    5   9       2.1
908  117  BH  15.854418   D229C  2023-04-24    D2  29    C    4  13       8.1
909  117  BH  15.854418   D229D  2023-04-24    D2  29    D    2   9      10.9
910  117  BH  15.854418   D235D  2023-04-24    D2  35    D    3   6      16.3
911  117  BH  15.854418   F115A  2023-04-24    F1  15    A    6   9      33.1
912  117  BH  15.854418   F229B  2023-04-24    F2  29    B    4   5       9.1
913  117  BH  15.854418   F240B  2023-04-24    F2  40    B    5  12      10.1
914  117  BH  15.854418   F235D  2023-04-24    F2  35    D    2   1       1.8
915  117  BH  15.854418    H14B  2023-04-24    H1   4    B    2  12      14.9
916  117  BH  15.854418    H18B  2023-04-24    H1   8    B    3   9      21.4
917  117  BH  15.854418   H114A  2023-04-24    H1  14    A    5   3      17.5
918  117  BH  15.854418   H224B  2023-04-24    H2  24    B    5  29      28.3
919  117  BH  15.854418   H233B  2023-04-24    H2  33    B    4   7      13.7
920  117  BH  15.854418   H110C  2023-04-24    H1  10    C    6  11      24.4
921  117  BH  15.854418   H223C  2023-04-24    H2  23    C    7  22      12.2
922  117  BH  15.854418    J14B  2023-04-24    J1   4    B    2   4       5.2
923  117  BH  15.854418   J119B  2023-04-24    J1  19    B    3  11      11.9
924  117  BH  15.854418   J233A  2023-04-24    J2  33    A    3   2       7.9
925  117  BH  15.854418   J237B  2023-04-24    J2  37    B    1   7      25.3
926  117  BH  15.854418   J110C  2023-04-24    J1  10    C    3   1      14.8
927  117  BH  15.854418   J119C  2023-04-24    J1  19    C    5  20      28.3
928  117  BH  15.854418   J223D  2023-04-24    J2  23    D    4   9       8.4
929  117  BH  15.854418   J225D  2023-04-24    J2  25    D    2   5      14.1
930  117  BH  15.854418   J230C  2023-04-24    J2  30    C    4  14      15.3
931  117  BH  15.854418   J231D  2023-04-24    J2  31    D    7  14      10.6
932  117  BH  15.854418   J232D  2023-04-24    J2  32    D    3  12      19.2
933  117  BH  15.854418   L110A  2023-04-24    L1  10    A    2  14      14.4
934  117  BH  15.854418    L14D  2023-04-24    L1   4    D    4  10      20.6
935  117  BH  15.854418   L113C  2023-04-24    L1  13    C    3  13      27.7
936  117  BH  15.854418   L113D  2023-04-24    L1  13    D    5  22      17.5
937  117  BH  15.854418   L228D  2023-04-24    L2  28    D    2   7       6.9
938  117  BH  15.854418   L229D  2023-04-24    L2  29    D    3  14      11.8
939  117  BH  15.854418   L240D  2023-04-24    L2  40    D    7  13      14.1
940  117  BH  15.854418   L241D  2023-04-24    L2  41    D    4   1      12.3
941  117 DPR  25.628088    F16C  2023-04-24    F1   6    C    6   9      16.0
942  117 DPR  25.628088   H117B  2023-04-24    H1  17    B    5   2      42.3
943  117 DPR  25.628088   H225C  2023-04-24    H2  25    C    2   6      36.1
944  117 DPR  25.628088    J16B  2023-04-24    J1   6    B    6  10      26.7
945  118  CC  18.755493    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
946  118  BH  16.905609    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
947  118 DPR  26.637579    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
948  119  CC  19.343950    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
949  119  BH  17.323663    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
950  119 DPR  27.839494    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
951  120  CC  19.966314    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
952  120  BH  17.810464    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
953  120 DPR  28.684801    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
954  121  CC  20.728691    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
955  121  BH  18.033895    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
956  121 DPR  29.772221    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
957  122  CC  21.559491    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
958  122  BH  19.186089    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
959  122 DPR  30.705082    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
960  123  CC  22.106187    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
961  123  BH  19.234210    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
962  123 DPR  32.116078    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
963  124  CC  23.186015   H111B  2023-05-01    H1  11    B    7   2      15.5
964  124  CC  23.186015   H115B  2023-05-01    H1  15    B    3   3      34.7
965  124  CC  23.186015   H116A  2023-05-01    H1  16    A    4  10      13.5
966  124  BH  20.139155   D226B  2023-05-01    D2  26    B    5  15      20.5
967  124  BH  20.139155   D230B  2023-05-01    D2  30    B    4   3      17.2
968  124  BH  20.139155   D119D  2023-05-01    D1  19    D    5   9       3.1
969  124  BH  20.139155   D229C  2023-05-01    D2  29    C    4  13      10.8
970  124  BH  20.139155   D229D  2023-05-01    D2  29    D    2   9      13.1
971  124  BH  20.139155   D235D  2023-05-01    D2  35    D    3   6       3.9
972  124  BH  20.139155   F115A  2023-05-01    F1  15    A    6   9      42.5
973  124  BH  20.139155   F229B  2023-05-01    F2  29    B    4   5      12.5
974  124  BH  20.139155   F240B  2023-05-01    F2  40    B    5  12      16.5
975  124  BH  20.139155   F235D  2023-05-01    F2  35    D    2   1      10.2
976  124  BH  20.139155    H14B  2023-05-01    H1   4    B    2  12      18.6
977  124  BH  20.139155    H18B  2023-05-01    H1   8    B    3   9      30.5
978  124  BH  20.139155   H114A  2023-05-01    H1  14    A    5   3      19.5
979  124  BH  20.139155   H224B  2023-05-01    H2  24    B    5  29      38.4
980  124  BH  20.139155   H233B  2023-05-01    H2  33    B    4   7      19.1
981  124  BH  20.139155   H110C  2023-05-01    H1  10    C    6  11      29.9
982  124  BH  20.139155   H223C  2023-05-01    H2  23    C    7  22      12.1
983  124  BH  20.139155    J14B  2023-05-01    J1   4    B    2   4      11.5
984  124  BH  20.139155   J119B  2023-05-01    J1  19    B    3  11      18.4
985  124  BH  20.139155   J233A  2023-05-01    J2  33    A    3   2       9.5
986  124  BH  20.139155   J237B  2023-05-01    J2  37    B    1   7      33.0
987  124  BH  20.139155   J110C  2023-05-01    J1  10    C    3   1      31.5
988  124  BH  20.139155   J119C  2023-05-01    J1  19    C    5  20      44.3
989  124  BH  20.139155   J223D  2023-05-01    J2  23    D    4   9      10.2
990  124  BH  20.139155   J225D  2023-05-01    J2  25    D    2   5      17.0
991  124  BH  20.139155   J230C  2023-05-01    J2  30    C    4  14      23.5
992  124  BH  20.139155   J231D  2023-05-01    J2  31    D    7  14      10.5
993  124  BH  20.139155   J232D  2023-05-01    J2  32    D    3  12      23.5
994  124  BH  20.139155   L110A  2023-05-01    L1  10    A    2  14      21.0
995  124  BH  20.139155    L14D  2023-05-01    L1   4    D    4  10      24.7
996  124  BH  20.139155   L113C  2023-05-01    L1  13    C    3  13      50.6
997  124  BH  20.139155   L113D  2023-05-01    L1  13    D    5  22      21.1
998  124  BH  20.139155   L228D  2023-05-01    L2  28    D    2   7       7.5
999  124  BH  20.139155   L229D  2023-05-01    L2  29    D    3  14      16.5
1000 124  BH  20.139155   L240D  2023-05-01    L2  40    D    7  13      15.0
1001 124  BH  20.139155   L241D  2023-05-01    L2  41    D    4   1      16.5
1002 124 DPR  33.312686    F16C  2023-05-01    F1   6    C    6   9      17.0
1003 124 DPR  33.312686   H117B  2023-05-01    H1  17    B    5   2      42.4
1004 124 DPR  33.312686    J16B  2023-05-01    J1   6    B    6  10      23.5
1005 125  CC  24.206748    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1006 125  BH  20.873723    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1007 125 DPR  34.263355    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1008 126  CC  24.467252    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1009 126  BH  21.764547    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1010 126 DPR  35.509791    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1011 127  CC  25.460823    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1012 127  BH  22.319869    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1013 127 DPR  36.888536    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1014 128  CC  26.273912    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1015 128  BH  23.336936    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1016 128 DPR  38.081594    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1017 129  CC  26.946204    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1018 129  BH  23.642292    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1019 129 DPR  39.689037    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1020 130  CC  27.719773    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1021 130  BH  24.694624    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1022 130 DPR  40.871805    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1023 131  CC  29.295713   H116A  2023-05-08    H1  16    A    4  10      15.1
1024 131  BH  25.454814   D226B  2023-05-08    D2  26    B    5  15      23.1
1025 131  BH  25.454814   D230B  2023-05-08    D2  30    B    4   3      22.1
1026 131  BH  25.454814   D119D  2023-05-08    D1  19    D    5   9       3.1
1027 131  BH  25.454814   D229C  2023-05-08    D2  29    C    4  13      13.3
1028 131  BH  25.454814   D229D  2023-05-08    D2  29    D    2   9      15.8
1029 131  BH  25.454814   D235D  2023-05-08    D2  35    D    3   6      35.2
1030 131  BH  25.454814   F115A  2023-05-08    F1  15    A    6   9      51.3
1031 131  BH  25.454814   F229B  2023-05-08    F2  29    B    4   5      14.6
1032 131  BH  25.454814   F240B  2023-05-08    F2  40    B    5  12      18.7
1033 131  BH  25.454814   F235D  2023-05-08    F2  35    D    2   1      12.2
1034 131  BH  25.454814    H14B  2023-05-08    H1   4    B    2  12      20.8
1035 131  BH  25.454814    H18B  2023-05-08    H1   8    B    3   9      36.1
1036 131  BH  25.454814   H114A  2023-05-08    H1  14    A    5   3      21.6
1037 131  BH  25.454814   H224B  2023-05-08    H2  24    B    5  29      48.3
1038 131  BH  25.454814   H233B  2023-05-08    H2  33    B    4   7      20.8
1039 131  BH  25.454814   H110C  2023-05-08    H1  10    C    6  11      34.0
1040 131  BH  25.454814   H223C  2023-05-08    H2  23    C    7  22      12.1
1041 131  BH  25.454814    J14B  2023-05-08    J1   4    B    2   4       4.4
1042 131  BH  25.454814   J119B  2023-05-08    J1  19    B    3  11      18.6
1043 131  BH  25.454814   J233A  2023-05-08    J2  33    A    3   2       5.7
1044 131  BH  25.454814   J237B  2023-05-08    J2  37    B    1   7      39.2
1045 131  BH  25.454814   J110C  2023-05-08    J1  10    C    3   1      33.1
1046 131  BH  25.454814   J119C  2023-05-08    J1  19    C    5  20      60.5
1047 131  BH  25.454814   J223D  2023-05-08    J2  23    D    4   9      13.2
1048 131  BH  25.454814   J225D  2023-05-08    J2  25    D    2   5      20.4
1049 131  BH  25.454814   J230C  2023-05-08    J2  30    C    4  14      24.6
1050 131  BH  25.454814   J231D  2023-05-08    J2  31    D    7  14       8.9
1051 131  BH  25.454814   J232D  2023-05-08    J2  32    D    3  12      24.0
1052 131  BH  25.454814   L110A  2023-05-08    L1  10    A    2  14      23.3
1053 131  BH  25.454814    L14D  2023-05-08    L1   4    D    4  10      27.3
1054 131  BH  25.454814   L113C  2023-05-08    L1  13    C    3  13      65.9
1055 131  BH  25.454814   L113D  2023-05-08    L1  13    D    5  22      21.2
1056 131  BH  25.454814   L228D  2023-05-08    L2  28    D    2   7       5.0
1057 131  BH  25.454814   L229D  2023-05-08    L2  29    D    3  14      17.8
1058 131  BH  25.454814   L240D  2023-05-08    L2  40    D    7  13      14.1
1059 131  BH  25.454814   L241D  2023-05-08    L2  41    D    4   1      11.5
1060 131 DPR  42.610789   H117B  2023-05-08    H1  17    B    5   2      42.2
1061 132  CC  30.199515    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1062 132  BH  26.121498    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1063 132 DPR  43.627068    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1064 133  CC  30.892772    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1065 133  BH  27.009168    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1066 133 DPR  45.399532    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1067 134  CC  31.795332    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1068 134  BH  28.020810    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1069 134 DPR  46.776527    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1070 135  CC  32.843026    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1071 135  BH  28.529934    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1072 135 DPR  48.191192    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1073 136  CC  33.815855    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1074 136  BH  29.214217    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1075 136 DPR  49.759878    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1076 137  CC  35.093420    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1077 137  BH  30.225219    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1078 137 DPR  51.629990    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1079 138  CC  35.718807    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1080 138  BH  31.528001    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1081 138 DPR  52.758754    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1082 139  CC  36.655352    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1083 139  BH  32.586632    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1084 139 DPR  54.413747    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1085 140  CC  37.835709    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1086 140  BH  32.735022    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1087 140 DPR  56.268015    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1088 141  CC  38.743501    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1089 141  BH  33.744373    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1090 141 DPR  57.678101    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1091 142  CC  40.341936    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1092 142  BH  35.199717    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1093 142 DPR  58.914588    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1094 143  CC  40.991558    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1095 143  BH  35.402047    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1096 143 DPR  60.672022    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1097 144  CC  42.011889    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1098 144  BH  36.684959    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1099 144 DPR  62.590425    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1100 145  CC  43.428819    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1101 145  BH  37.452258    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1102 145 DPR  64.170865    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1103 146  CC  44.483971    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1104 146  BH  38.749361    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1105 146 DPR  66.013778    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1106 147  CC  45.344375    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1107 147  BH  39.267125    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1108 147 DPR  67.292749    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1109 148  CC  46.635660    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1110 148  BH  40.325140    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1111 148 DPR  69.129043    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1112 149  CC  47.444503    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1113 149  BH  40.857720    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1114 149 DPR  70.870036    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1115 150  CC  48.758416    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1116 150  BH  41.951961    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1117 150 DPR  72.417466    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1118 151  CC  49.318200    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1119 151  BH  43.201930    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1120 151 DPR  73.797206    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1121 152  CC  50.829998    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1122 152  BH  43.990472    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1123 152 DPR  75.570837    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1124 153  CC  51.754827    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1125 153  BH  44.748020    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1126 153 DPR  77.512828    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1127 154  CC  52.533355    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1128 154  BH  45.424700    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1129 154 DPR  78.861017    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1130 155  CC  53.760669    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1131 155  BH  46.383377    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1132 155 DPR  80.207348    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1133 156  CC  54.646302    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1134 156  BH  47.357999    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1135 156 DPR  82.094322    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1136 157  CC  55.778727    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1137 157  BH  48.310541    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1138 157 DPR  83.272250    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1139 158  CC  56.713847    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1140 158  BH  49.135436    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1141 158 DPR  84.788434    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1142 159  CC  57.612095    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1143 159  BH  49.734185   D226B  2023-06-05    D2  26    B    5  15      77.4
1144 159  BH  49.734185   D230B  2023-06-05    D2  30    B    4   3      52.1
1145 159  BH  49.734185   D229C  2023-06-05    D2  29    C    4  13      23.7
1146 159  BH  49.734185   D229D  2023-06-05    D2  29    D    2   9      34.5
1147 159  BH  49.734185   D235D  2023-06-05    D2  35    D    3   6      80.5
1148 159  BH  49.734185   F115A  2023-06-05    F1  15    A    6   9      88.9
1149 159  BH  49.734185   F229B  2023-06-05    F2  29    B    4   5      26.3
1150 159  BH  49.734185   F240B  2023-06-05    F2  40    B    5  12      47.0
1151 159  BH  49.734185   F235D  2023-06-05    F2  35    D    2   1      37.5
1152 159  BH  49.734185    H14B  2023-06-05    H1   4    B    2  12      28.2
1153 159  BH  49.734185    H18B  2023-06-05    H1   8    B    3   9      76.1
1154 159  BH  49.734185   H114A  2023-06-05    H1  14    A    5   3      19.9
1155 159  BH  49.734185   H224B  2023-06-05    H2  24    B    5  29      70.8
1156 159  BH  49.734185   H233B  2023-06-05    H2  33    B    4   7      21.1
1157 159  BH  49.734185   H110C  2023-06-05    H1  10    C    6  11      44.8
1158 159  BH  49.734185   J119B  2023-06-05    J1  19    B    3  11      31.9
1159 159  BH  49.734185   J237B  2023-06-05    J2  37    B    1   7      39.9
1160 159  BH  49.734185   J110C  2023-06-05    J1  10    C    3   1      70.0
1161 159  BH  49.734185   J119C  2023-06-05    J1  19    C    5  20     102.4
1162 159  BH  49.734185   J223D  2023-06-05    J2  23    D    4   9      51.4
1163 159  BH  49.734185   J225D  2023-06-05    J2  25    D    2   5      71.0
1164 159  BH  49.734185   J230C  2023-06-05    J2  30    C    4  14      32.3
1165 159  BH  49.734185   J232D  2023-06-05    J2  32    D    3  12      27.8
1166 159  BH  49.734185   L110A  2023-06-05    L1  10    A    2  14      61.7
1167 159  BH  49.734185    L14D  2023-06-05    L1   4    D    4  10      77.6
1168 159  BH  49.734185   L113C  2023-06-05    L1  13    C    3  13     120.6
1169 159  BH  49.734185   L113D  2023-06-05    L1  13    D    5  22      58.9
1170 159  BH  49.734185   L229D  2023-06-05    L2  29    D    3  14      31.4
1171 159  BH  49.734185   L240D  2023-06-05    L2  40    D    7  13      19.7
1172 159  BH  49.734185   L241D  2023-06-05    L2  41    D    4   1      49.4
1173 159 DPR  86.471080    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1174 160  CC  58.997514    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1175 160  BH  50.688490    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1176 160 DPR  87.655038    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1177 161  CC  59.689022    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1178 161  BH  51.206258    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1179 161 DPR  89.275447    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1180 162  CC  60.444328    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1181 162  BH  52.396706    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1182 162 DPR  90.660921    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1183 163  CC  61.160444    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1184 163  BH  53.058676    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1185 163 DPR  91.815749    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1186 164  CC  62.019354    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1187 164  BH  53.299503    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1188 164 DPR  93.246374    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1189 165  CC  63.336032    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1190 165  BH  54.077978    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1191 165 DPR  94.467308    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1192 166  CC  63.671354    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1193 166  BH  54.991396    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1194 166 DPR  95.394941    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1195 167  CC  64.292896    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1196 167  BH  55.804415    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1197 167 DPR  97.088646    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1198 168  CC  65.367758    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1199 168  BH  56.450933    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1200 168 DPR  97.962900    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1201 169  CC  66.199357    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1202 169  BH  56.461130    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1203 169 DPR  99.241612    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1204 170  CC  66.700881    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1205 170  BH  57.656488    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1206 170 DPR  99.902797    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1207 171  CC  67.068760    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1208 171  BH  58.252088    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1209 171 DPR 100.985188    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1210 172  CC  67.825510    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1211 172  BH  58.319256    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1212 172 DPR 101.710708    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1213 173  CC  68.525348    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1214 173  BH  59.023606    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1215 173 DPR 103.202053    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1216 174  CC  69.229135    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1217 174  BH  59.550697    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1218 174 DPR 103.745338    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1219 175  CC  69.683937    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1220 175  BH  60.219075    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1221 175 DPR 104.540294    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1222 176  CC  69.921073    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1223 176  BH  60.619012    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1224 176 DPR 105.111218    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1225 177  CC  70.664155    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1226 177  BH  60.887290    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1227 177 DPR 106.243118    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1228 178  CC  70.999329    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1229 178  BH  61.350592    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1230 178 DPR 106.743919    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1231 179  CC  71.537636    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1232 179  BH  61.359531    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1233 179 DPR 107.525110    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1234 180  CC  71.401900    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1235 180  BH  62.139732    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1236 180 DPR 108.406295    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1237 181  CC  72.118821    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1238 181  BH  62.339211    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1239 181 DPR 108.737296    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1240 182  CC  72.432180    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1241 182  BH  62.439519    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1242 182 DPR 108.999232    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1243 183  CC  73.088647    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1244 183  BH  62.426665    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1245 183 DPR 109.513524    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1246 184  CC  73.192719    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1247 184  BH  63.256294    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1248 184 DPR 110.091999    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1249 185  CC  73.545953    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1250 185  BH  63.411164    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1251 185 DPR 110.941857    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1252 186  CC  73.587069    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1253 186  BH  63.508790    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1254 186 DPR 111.006819    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1255 187  CC  74.275295    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1256 187  BH  63.817876   D226B  2023-07-03    D2  26    B    5  15      67.4
1257 187  BH  63.817876   D230B  2023-07-03    D2  30    B    4   3      76.3
1258 187  BH  63.817876   D229C  2023-07-03    D2  29    C    4  13      26.9
1259 187  BH  63.817876   D229D  2023-07-03    D2  29    D    2   9      42.7
1260 187  BH  63.817876   D235D  2023-07-03    D2  35    D    3   6      83.2
1261 187  BH  63.817876   F115A  2023-07-03    F1  15    A    6   9     110.4
1262 187  BH  63.817876   F229B  2023-07-03    F2  29    B    4   5      34.3
1263 187  BH  63.817876   F240B  2023-07-03    F2  40    B    5  12      70.9
1264 187  BH  63.817876   F235D  2023-07-03    F2  35    D    2   1      67.8
1265 187  BH  63.817876    H18B  2023-07-03    H1   8    B    3   9      81.6
1266 187  BH  63.817876   H114A  2023-07-03    H1  14    A    5   3      18.7
1267 187  BH  63.817876   H224B  2023-07-03    H2  24    B    5  29      69.0
1268 187  BH  63.817876   H233B  2023-07-03    H2  33    B    4   7      18.0
1269 187  BH  63.817876   H110C  2023-07-03    H1  10    C    6  11      45.1
1270 187  BH  63.817876   J119B  2023-07-03    J1  19    B    3  11      32.1
1271 187  BH  63.817876   J110C  2023-07-03    J1  10    C    3   1      67.0
1272 187  BH  63.817876   J119C  2023-07-03    J1  19    C    5  20     101.0
1273 187  BH  63.817876   J223D  2023-07-03    J2  23    D    4   9      49.0
1274 187  BH  63.817876   J225D  2023-07-03    J2  25    D    2   5      82.0
1275 187  BH  63.817876   J230C  2023-07-03    J2  30    C    4  14      35.0
1276 187  BH  63.817876   J232D  2023-07-03    J2  32    D    3  12      29.0
1277 187  BH  63.817876   L110A  2023-07-03    L1  10    A    2  14      89.5
1278 187  BH  63.817876    L14D  2023-07-03    L1   4    D    4  10      92.5
1279 187  BH  63.817876   L113C  2023-07-03    L1  13    C    3  13      82.0
1280 187  BH  63.817876   L113D  2023-07-03    L1  13    D    5  22     113.0
1281 187  BH  63.817876   L229D  2023-07-03    L2  29    D    3  14      28.0
1282 187  BH  63.817876   L240D  2023-07-03    L2  40    D    7  13      17.0
1283 187  BH  63.817876   L241D  2023-07-03    L2  41    D    4   1      67.5
1284 187 DPR 111.421777    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1285 188  CC  74.025731    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1286 188  BH  63.524032    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1287 188 DPR 111.581988    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1288 189  CC  74.212877    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1289 189  BH  63.809744    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1290 189 DPR 111.752591    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1291 190  CC  74.581388    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1292 190  BH  64.493028    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1293 190 DPR 112.219144    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1294 191  CC  74.711978    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1295 191  BH  64.605901    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1296 191 DPR 112.448065    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1297 192  CC  74.716042    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1298 192  BH  64.506922    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1299 192 DPR 112.976009    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1300 193  CC  74.993750    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1301 193  BH  64.501967    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1302 193 DPR 113.012222    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1303 194  CC  74.777051    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1304 194  BH  64.683957   D226B  2023-07-10    D2  26    B    5  15      70.0
1305 194  BH  64.683957   D230B  2023-07-10    D2  30    B    4   3      73.0
1306 194  BH  64.683957   D229C  2023-07-10    D2  29    C    4  13      26.5
1307 194  BH  64.683957   D229D  2023-07-10    D2  29    D    2   9      43.0
1308 194  BH  64.683957   D235D  2023-07-10    D2  35    D    3   6      85.8
1309 194  BH  64.683957   F115A  2023-07-10    F1  15    A    6   9     113.2
1310 194  BH  64.683957   F229B  2023-07-10    F2  29    B    4   5      33.2
1311 194  BH  64.683957   F240B  2023-07-10    F2  40    B    5  12      71.6
1312 194  BH  64.683957   F235D  2023-07-10    F2  35    D    2   1      68.1
1313 194  BH  64.683957    H18B  2023-07-10    H1   8    B    3   9      82.2
1314 194  BH  64.683957   H224B  2023-07-10    H2  24    B    5  29      68.1
1315 194  BH  64.683957   H233B  2023-07-10    H2  33    B    4   7      23.5
1316 194  BH  64.683957   J223D  2023-07-10    J2  23    D    4   9      49.6
1317 194  BH  64.683957   J225D  2023-07-10    J2  25    D    2   5      77.0
1318 194  BH  64.683957   J230C  2023-07-10    J2  30    C    4  14      35.1
1319 194  BH  64.683957   J232D  2023-07-10    J2  32    D    3  12      30.5
1320 194  BH  64.683957   L110A  2023-07-10    L1  10    A    2  14      89.4
1321 194  BH  64.683957    L14D  2023-07-10    L1   4    D    4  10      95.6
1322 194  BH  64.683957   L113C  2023-07-10    L1  13    C    3  13     130.0
1323 194  BH  64.683957   L113D  2023-07-10    L1  13    D    5  22      96.1
1324 194  BH  64.683957   L240D  2023-07-10    L2  40    D    7  13      20.1
1325 194  BH  64.683957   L241D  2023-07-10    L2  41    D    4   1      68.4
1326 194 DPR 113.161107    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1327 195  CC  75.346477    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1328 195  BH  64.824782    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1329 195 DPR 113.571667    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1330 196  CC  75.690129    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1331 196  BH  64.184288    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1332 196 DPR 113.636502    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1333 197  CC  75.238847    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1334 197  BH  64.968679    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1335 197 DPR 113.701174    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1336 198  CC  75.632811    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1337 198  BH  65.323160    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1338 198 DPR 113.730190    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1339 199  CC  75.440489    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1340 199  BH  65.135340    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1341 199 DPR 113.877032    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1342 200  CC  75.936859    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1343 200  BH  65.293806    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1344 200 DPR 114.136976    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1345 201  CC  75.790892    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1346 201  BH  65.248888   D226B  2023-07-17    D2  26    B    5  15      68.2
1347 201  BH  65.248888   D230B  2023-07-17    D2  30    B    4   3      77.8
1348 201  BH  65.248888   D229C  2023-07-17    D2  29    C    4  13      25.7
1349 201  BH  65.248888   D229D  2023-07-17    D2  29    D    2   9      44.2
1350 201  BH  65.248888   D235D  2023-07-17    D2  35    D    3   6      83.6
1351 201  BH  65.248888   F115A  2023-07-17    F1  15    A    6   9     113.2
1352 201  BH  65.248888   F229B  2023-07-17    F2  29    B    4   5      34.3
1353 201  BH  65.248888   F240B  2023-07-17    F2  40    B    5  12      76.3
1354 201  BH  65.248888   F235D  2023-07-17    F2  35    D    2   1      68.8
1355 201  BH  65.248888    H18B  2023-07-17    H1   8    B    3   9      83.7
1356 201  BH  65.248888   H224B  2023-07-17    H2  24    B    5  29      70.3
1357 201  BH  65.248888   H233B  2023-07-17    H2  33    B    4   7      22.9
1358 201  BH  65.248888   J223D  2023-07-17    J2  23    D    4   9      50.8
1359 201  BH  65.248888   J225D  2023-07-17    J2  25    D    2   5      83.5
1360 201  BH  65.248888   J230C  2023-07-17    J2  30    C    4  14      41.6
1361 201  BH  65.248888   J232D  2023-07-17    J2  32    D    3  12      30.7
1362 201  BH  65.248888   L110A  2023-07-17    L1  10    A    2  14      90.0
1363 201  BH  65.248888    L14D  2023-07-17    L1   4    D    4  10      73.2
1364 201  BH  65.248888   L113C  2023-07-17    L1  13    C    3  13     107.6
1365 201  BH  65.248888   L241D  2023-07-17    L2  41    D    4   1      37.4
1366 201 DPR 114.148383    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1367 202  CC  75.600649    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1368 202  BH  65.001724    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1369 202 DPR 114.008513    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1370 203  CC  75.738017    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1371 203  BH  65.327347    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1372 203 DPR 113.992887    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1373 204  CC  75.757523    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1374 204  BH  65.443451    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1375 204 DPR 114.438949    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1376 205  CC  76.205540    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1377 205  BH  65.590407    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1378 205 DPR 114.387549    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1379 206  CC  76.073916    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1380 206  BH  65.498734    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1381 206 DPR 114.400217    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1382 207  CC  76.154203    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1383 207  BH  65.295931    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1384 207 DPR 114.260722    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1385 208  CC  76.149654    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1386 208  BH  65.426895   D226B  2023-07-24    D2  26    B    5  15      68.8
1387 208  BH  65.426895   D230B  2023-07-24    D2  30    B    4   3      79.3
1388 208  BH  65.426895   D229C  2023-07-24    D2  29    C    4  13      26.4
1389 208  BH  65.426895   D229D  2023-07-24    D2  29    D    2   9      44.8
1390 208  BH  65.426895   D235D  2023-07-24    D2  35    D    3   6      86.0
1391 208  BH  65.426895   F115A  2023-07-24    F1  15    A    6   9     109.9
1392 208  BH  65.426895   F229B  2023-07-24    F2  29    B    4   5      36.2
1393 208  BH  65.426895   F240B  2023-07-24    F2  40    B    5  12      75.1
1394 208  BH  65.426895   F235D  2023-07-24    F2  35    D    2   1      69.3
1395 208  BH  65.426895    H18B  2023-07-24    H1   8    B    3   9      81.7
1396 208  BH  65.426895   H224B  2023-07-24    H2  24    B    5  29      68.5
1397 208  BH  65.426895   H233B  2023-07-24    H2  33    B    4   7      22.8
1398 208  BH  65.426895   J225D  2023-07-24    J2  25    D    2   5      81.5
1399 208  BH  65.426895   J230C  2023-07-24    J2  30    C    4  14      32.8
1400 208  BH  65.426895   J232D  2023-07-24    J2  32    D    3  12      30.2
1401 208  BH  65.426895   L110A  2023-07-24    L1  10    A    2  14      91.2
1402 208  BH  65.426895    L14D  2023-07-24    L1   4    D    4  10      93.3
1403 208  BH  65.426895   L113C  2023-07-24    L1  13    C    3  13      89.8
1404 208  BH  65.426895   L113D  2023-07-24    L1  13    D    5  22     112.1
1405 208  BH  65.426895   L241D  2023-07-24    L2  41    D    4   1      69.6
1406 208 DPR 114.288121    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
     herbivory     pred2
1         <NA>  4.148911
2         <NA>  3.842604
3         <NA>  3.993838
4            Y  4.105850
5         <NA>  4.148758
6            Y  3.972188
7            Y  4.243553
8         <NA>  4.554645
9         <NA>  4.167385
10        <NA>  3.979970
11           Y  4.095179
12        <NA>  4.029672
13           Y  4.077110
14        <NA>  4.099911
15        <NA>  4.014169
16        <NA>  3.940888
17           Y  3.655683
18        <NA>  3.936292
19           Y  3.696057
20           Y  3.947541
21        <NA>  4.109962
22        <NA>  4.197873
23        <NA>  3.945455
24        <NA>  3.681728
25        <NA>  4.100692
26        <NA>  4.037502
27        <NA>  4.196459
28           Y  3.734466
29        <NA>  4.073610
30        <NA>  4.001050
31           Y  4.277402
32        <NA>  3.911090
33        <NA>  4.278304
34        <NA>  3.735580
35        <NA>  4.073783
36           Y  4.222813
37        <NA>  3.823875
38        <NA>  4.044276
39        <NA>  3.737831
40        <NA>  3.815263
41        <NA>  3.657853
42        <NA>  4.217353
43        <NA>  4.022127
44        <NA>  4.075006
45        <NA>  3.960845
46        <NA>  4.108784
47        <NA>  4.415648
48        <NA>  4.034198
49        <NA>  4.037007
50        <NA>  4.177130
51           Y  3.722831
52        <NA>  3.814003
53        <NA>  3.802641
54        <NA>  3.990997
55        <NA>  4.090552
56           Y  4.313665
57        <NA>  3.984565
58        <NA>  4.173173
59        <NA>  3.920891
60        <NA>  4.113407
61           Y  4.175665
62        <NA>  4.338456
63        <NA>  3.594304
64        <NA>  3.892561
65        <NA>  3.820037
66        <NA>  4.097466
67        <NA>  4.138569
68        <NA>  4.263391
69        <NA>  3.779162
70           Y  3.895233
71        <NA>  3.979383
72        <NA>  4.350284
73           Y  3.915123
74        <NA>  3.836081
75           Y  3.982300
76           Y  4.092511
77           Y  4.149602
78           Y  4.002127
79        <NA>  3.746798
80           Y  3.957699
81        <NA>  3.951152
82        <NA>  3.925522
83        <NA>  4.120181
84           Y  3.865639
85           Y  3.865019
86           Y  4.043930
87        <NA>  4.153882
88        <NA>  4.037575
89           Y  3.852776
90        <NA>  3.965145
91           Y  3.835227
92           Y  3.914002
93           Y  4.081494
94           Y  4.210114
95        <NA>  4.189636
96           Y  4.515329
97           Y  3.685171
98        <NA>  4.285534
99        <NA>  3.870562
100          Y  4.075742
101       <NA>  3.614083
102          Y  4.345371
103       <NA>  3.960849
104       <NA>  4.195113
105       <NA>  4.039424
106          Y  4.219755
107       <NA>  4.335103
108       <NA>  4.005785
109          Y  3.970489
110       <NA>  4.206630
111       <NA>  3.926732
112       <NA>  3.914115
113       <NA>  3.911567
114          Y  4.002342
115          Y  3.599615
116       <NA>  3.843831
117       <NA>  4.044912
118       <NA>  3.958283
119          Y  3.618131
120       <NA>  4.130432
121          Y  3.874597
122       <NA>  4.141417
123       <NA>  4.330390
124       <NA>  4.141104
125       <NA>  4.315238
126          Y  4.077982
127          Y  3.866668
128          Y  3.981769
129       <NA>  3.866832
130          Y  3.967090
131       <NA>  3.940313
132          Y  3.902978
133          Y  3.952778
134       <NA>  4.285802
135          Y  3.953461
136       <NA>  4.093160
137          Y  4.194416
138          Y  4.119049
139       <NA>  4.043877
140       <NA>  4.243309
141          Y  3.802787
142       <NA>  4.298092
143       <NA>  3.965468
144       <NA>  4.046314
145          Y  4.005669
146       <NA>  3.747772
147          Y  4.155780
148       <NA>  3.675396
149       <NA>  4.253410
150       <NA>  4.071503
151       <NA>  3.812498
152       <NA>  4.078447
153       <NA>  4.136399
154       <NA>  4.275909
155          Y  4.434647
156       <NA>        NA
157       <NA>        NA
158       <NA>        NA
159       <NA>        NA
160       <NA>        NA
161       <NA>        NA
162       <NA>        NA
163       <NA>        NA
164       <NA>        NA
165       <NA>        NA
166       <NA>        NA
167       <NA>        NA
168       <NA>        NA
169       <NA>        NA
170       <NA>        NA
171       <NA>        NA
172       <NA>        NA
173       <NA>        NA
174       <NA>        NA
175       <NA>        NA
176       <NA>        NA
177       <NA>        NA
178       <NA>        NA
179       <NA>        NA
180       <NA>        NA
181       <NA>        NA
182       <NA>        NA
183       <NA>        NA
184       <NA>        NA
185       <NA>        NA
186       <NA>        NA
187       <NA>        NA
188       <NA>        NA
189       <NA>        NA
190       <NA>        NA
191       <NA>        NA
192       <NA>        NA
193       <NA>        NA
194       <NA>        NA
195       <NA>        NA
196       <NA>        NA
197       <NA>        NA
198       <NA>        NA
199       <NA>        NA
200       <NA>        NA
201       <NA>        NA
202       <NA>        NA
203       <NA>        NA
204       <NA>        NA
205       <NA>        NA
206       <NA>        NA
207       <NA>        NA
208       <NA>        NA
209       <NA>        NA
210       <NA>        NA
211       <NA>        NA
212       <NA>        NA
213       <NA>        NA
214       <NA>        NA
215       <NA>        NA
216       <NA>  4.048867
217       <NA>  4.177960
218       <NA>  4.469056
219       <NA>  4.168844
220       <NA>  3.723589
221       <NA>  4.330865
222       <NA>  4.090079
223       <NA>  3.487898
224       <NA>  4.180318
225       <NA>  4.218069
226       <NA>  4.118290
227       <NA>  3.798805
228       <NA>  4.112770
229       <NA>  4.246015
230       <NA>  4.365325
231       <NA>  3.992030
232       <NA>  4.386487
233          N  4.147818
234       <NA>  3.837858
235       <NA>  4.116936
236       <NA>  4.244860
237       <NA>  4.082601
238       <NA>  4.186298
239       <NA>  4.114825
240       <NA>  4.180244
241       <NA>  4.236157
242       <NA>  4.149765
243       <NA>  4.198047
244       <NA>  4.096834
245       <NA>  3.721473
246       <NA>  3.947717
247       <NA>  3.966014
248       <NA>  3.883421
249       <NA>  4.005281
250       <NA>  4.046464
251       <NA>  4.111995
252       <NA>  3.919265
253       <NA>  4.316375
254       <NA>  4.168178
255       <NA>  3.796882
256       <NA>  3.787917
257       <NA>  4.014536
258       <NA>  4.183102
259       <NA>  4.167752
260       <NA>  3.911518
261       <NA>  3.989009
262       <NA>  4.029773
263       <NA>  3.652474
264       <NA>  4.242688
265       <NA>  3.974901
266       <NA>  4.074099
267       <NA>  4.365486
268       <NA>  4.129809
269       <NA>  4.379698
270       <NA>  4.197630
271       <NA>  3.961687
272       <NA>  4.263864
273       <NA>  4.196400
274          N  4.155887
275       <NA>  4.145577
276       <NA>  3.961754
277       <NA>  3.780129
278       <NA>  3.980115
279       <NA>  3.914017
280       <NA>  4.188905
281       <NA>  4.320061
282       <NA>  3.958689
283       <NA>  4.129275
284       <NA>  4.053618
285       <NA>  4.131063
286       <NA>  4.379329
287       <NA>  3.858260
288       <NA>  4.166941
289       <NA>  4.081215
290       <NA>  3.852516
291       <NA>  4.211726
292       <NA>  4.296049
293       <NA>  3.801700
294       <NA>  4.283604
295       <NA>  4.201705
296       <NA>  4.070397
297       <NA>  4.074574
298       <NA>  4.179010
299       <NA>  4.010196
300       <NA>  4.051335
301       <NA>  4.059445
302       <NA>  4.290635
303       <NA>  3.842451
304       <NA>  3.943348
305       <NA>  3.845783
306       <NA>  4.274208
307       <NA>  3.831344
308       <NA>  4.153988
309       <NA>  4.295971
310       <NA>  4.088300
311       <NA>  4.334133
312       <NA>  4.631150
313       <NA>  4.158522
314       <NA>  3.970058
315       <NA>  4.032030
316       <NA>  4.281962
317       <NA>  4.136768
318       <NA>  4.337349
319       <NA>  4.027108
320       <NA>  3.997098
321       <NA>  4.039054
322       <NA>  4.262544
323       <NA>  3.945339
324       <NA>  4.484530
325       <NA>  4.173338
326       <NA>  3.965640
327       <NA>  3.988919
328       <NA>  4.009723
329       <NA>  4.476264
330       <NA>  4.313164
331       <NA>  3.905657
332       <NA>  4.480583
333       <NA>  4.275778
334       <NA>  4.156127
335       <NA>  3.902299
336       <NA>  4.211665
337       <NA>  4.249444
338       <NA>  4.002660
339       <NA>  4.443792
340       <NA>  4.009478
341       <NA>  4.462206
342          2  4.393135
343       <NA>        NA
344       <NA>        NA
345       <NA>        NA
346       <NA>        NA
347       <NA>        NA
348       <NA>        NA
349       <NA>        NA
350       <NA>        NA
351       <NA>        NA
352       <NA>        NA
353       <NA>        NA
354       <NA>        NA
355       <NA>        NA
356       <NA>        NA
357       <NA>        NA
358       <NA>        NA
359       <NA>        NA
360       <NA>        NA
361       <NA>        NA
362       <NA>        NA
363       <NA>        NA
364       <NA>        NA
365       <NA>        NA
366       <NA>        NA
367       <NA>        NA
368       <NA>        NA
369       <NA>        NA
370       <NA>        NA
371       <NA>        NA
372       <NA>        NA
373       <NA>        NA
374       <NA>        NA
375       <NA>        NA
376       <NA>        NA
377       <NA>        NA
378       <NA>        NA
379       <NA>        NA
380       <NA>        NA
381       <NA>        NA
382       <NA>  4.623939
383       <NA>  4.761528
384       <NA>  4.551041
385       <NA>  4.408479
386       <NA>  3.999645
387       <NA>  4.672601
388       <NA>  4.620095
389       <NA>  4.484841
390       <NA>  4.473938
391       <NA>  4.337483
392       <NA>  4.091913
393       <NA>  4.434084
394       <NA>  4.405738
395       <NA>  4.623842
396       <NA>  4.479289
397       <NA>  4.542694
398       <NA>  4.638358
399       <NA>  4.585283
400       <NA>  4.338908
401       <NA>  4.759371
402       <NA>  4.469668
403       <NA>  4.234865
404       <NA>  4.172339
405       <NA>  4.587297
406       <NA>  4.618472
407       <NA>  4.250877
408       <NA>  4.378571
409       <NA>  4.674958
410       <NA>  4.709782
411       <NA>  4.489206
412       <NA>  4.281229
413       <NA>  4.548151
414       <NA>  3.979702
415       <NA>  4.592712
416       <NA>  4.511183
417       <NA>  3.972744
418       <NA>  4.382486
419       <NA>  4.453472
420       <NA>  4.313156
421       <NA>  4.241630
422       <NA>  4.757375
423       <NA>  4.555152
424       <NA>  4.443488
425       <NA>  4.125467
426       <NA>  4.600950
427       <NA>  4.623972
428       <NA>  4.427216
429       <NA>  4.032581
430       <NA>  4.513716
431       <NA>  4.465592
432       <NA>  4.046029
433       <NA>  4.445083
434       <NA>  4.588436
435       <NA>  4.193445
436       <NA>  4.377026
437       <NA>  4.140356
438       <NA>  4.090222
439       <NA>  4.106166
440       <NA>  4.543175
441       <NA>  4.539635
442       <NA>  4.365450
443       <NA>  4.578484
444       <NA>  4.476942
445       <NA>  4.737652
446       <NA>  4.337525
447       <NA>  4.507555
448       <NA>  4.224565
449       <NA>  4.668668
450       <NA>  4.598945
451       <NA>  4.029952
452       <NA>  4.322487
453       <NA>  4.453681
454       <NA>  4.554868
455       <NA>  4.324717
456       <NA>  4.389765
457       <NA>  4.532365
458       <NA>  4.665198
459       <NA>  4.390023
460       <NA>  4.009658
461       <NA>  4.345473
462       <NA>  4.024504
463       <NA>  4.146161
464       <NA>  4.437604
465       <NA>  4.872679
466       <NA>  4.717533
467       <NA>  4.463176
468       <NA>  4.873331
469       <NA>  4.927516
470       <NA>  4.454619
471       <NA>  4.336828
472       <NA>  4.819753
473       <NA>  4.509873
474       <NA>  5.113462
475       <NA>  4.745904
476       <NA>        NA
477       <NA>        NA
478       <NA>        NA
479       <NA>        NA
480       <NA>        NA
481       <NA>        NA
482       <NA>        NA
483       <NA>        NA
484       <NA>        NA
485       <NA>        NA
486       <NA>        NA
487       <NA>        NA
488       <NA>        NA
489       <NA>        NA
490       <NA>        NA
491       <NA>        NA
492       <NA>        NA
493       <NA>        NA
494       <NA>        NA
495       <NA>        NA
496       <NA>        NA
497       <NA>        NA
498       <NA>        NA
499       <NA>        NA
500       <NA>        NA
501       <NA>        NA
502       <NA>        NA
503       <NA>        NA
504       <NA>        NA
505       <NA>        NA
506       <NA>        NA
507       <NA>        NA
508       <NA>        NA
509       <NA>        NA
510       <NA>        NA
511       <NA>        NA
512       <NA>        NA
513       <NA>        NA
514       <NA>        NA
515       <NA>  5.260101
516       <NA>  5.372108
517       <NA>  5.963634
518       <NA>  5.439538
519       <NA>  5.197355
520       <NA>  5.371875
521       <NA>  5.051516
522       <NA>  5.424802
523       <NA>  5.303648
524       <NA>  5.359812
525       <NA>  5.730514
526       <NA>  5.506673
527       <NA>  5.483910
528       <NA>  5.603551
529       <NA>  5.240583
530       <NA>  5.152138
531       <NA>  5.269413
532       <NA>  5.191016
533       <NA>  5.301069
534       <NA>  5.100961
535       <NA>  5.578705
536       <NA>  5.264534
537       <NA>  5.297179
538       <NA>  5.174086
539       <NA>  5.471559
540       <NA>  5.433759
541       <NA>  5.301555
542       <NA>  5.500230
543       <NA>  5.420361
544       <NA>  5.289902
545       <NA>  5.081176
546       <NA>  5.139862
547       <NA>  5.251753
548       <NA>  5.123939
549       <NA>  5.534674
550       <NA>  5.196862
551       <NA>  5.275141
552       <NA>  5.200976
553       <NA>  5.292952
554       <NA>  5.273062
555       <NA>  4.812851
556       <NA>  5.150580
557       <NA>  5.342203
558       <NA>  5.440135
559       <NA>  5.398504
560       <NA>  5.289795
561       <NA>  5.149202
562       <NA>  5.445818
563       <NA>  5.476427
564       <NA>  5.121850
565       <NA>  5.356923
566       <NA>  5.037571
567       <NA>  5.047232
568       <NA>  5.307702
569       <NA>  4.982859
570       <NA>  5.073155
571       <NA>  5.621993
572       <NA>  5.025220
573       <NA>  5.508355
574       <NA>  5.298746
575       <NA>  5.532122
576       <NA>  5.376421
577       <NA>  5.588459
578       <NA>  5.116145
579       <NA>  5.518743
580       <NA>  4.965381
581       <NA>  5.174467
582       <NA>  5.384502
583       <NA>  5.449907
584       <NA>  6.203346
585       <NA>  6.208705
586       <NA>  6.211303
587       <NA>  6.107454
588       <NA>  5.887029
589       <NA>  6.147627
590       <NA>  6.159257
591       <NA>  6.282652
592       <NA>  6.211095
593       <NA>        NA
594       <NA>        NA
595       <NA>        NA
596       <NA>        NA
597       <NA>        NA
598       <NA>        NA
599       <NA>        NA
600       <NA>        NA
601       <NA>        NA
602       <NA>        NA
603       <NA>        NA
604       <NA>        NA
605       <NA>        NA
606       <NA>        NA
607       <NA>        NA
608       <NA>        NA
609       <NA>        NA
610       <NA>        NA
611          N  6.342016
612       <NA>  6.367655
613       <NA>  6.591619
614       <NA>  6.454612
615       <NA>  6.772238
616          Y  6.395129
617       <NA>  6.668131
618       <NA>  6.600416
619          N  6.606587
620          N  6.296395
621       <NA>  6.605689
622       <NA>  6.176004
623       <NA>  6.248076
624       <NA>  6.080169
625       <NA>  6.182312
626       <NA>  6.061265
627       <NA>  5.927651
628       <NA>  5.938665
629          Y  6.079332
630       <NA>  5.936466
631       <NA>  6.317369
632       <NA>  6.393786
633       <NA>  6.017141
634       <NA>  6.079569
635       <NA>  6.171721
636       <NA>  6.510651
637       <NA>  6.161970
638       <NA>  6.515056
639       <NA>  5.913427
640          Y  6.046823
641          Y  6.079364
642       <NA>  5.898380
643          Y  5.839459
644       <NA>  6.089219
645       <NA>  6.236456
646       <NA>  5.875362
647       <NA>  6.213592
648       <NA>  6.131633
649       <NA>  5.943020
650       <NA>  6.163403
651       <NA>  5.734252
652       <NA>  6.485041
653       <NA>  6.158781
654       <NA>  6.416736
655       <NA>  5.804046
656       <NA>  5.964443
657       <NA>  5.981661
658       <NA>  5.833019
659       <NA>  6.149303
660       <NA>  6.388487
661       <NA>  6.159071
662       <NA>  5.923625
663       <NA>  6.269249
664          Y  6.018549
665       <NA>  6.024325
666       <NA>  5.801129
667       <NA>  7.918966
668       <NA>  7.846969
669       <NA>  7.482148
670       <NA>  7.578115
671       <NA>  7.737885
672       <NA>  7.526194
673       <NA>  7.353018
674          N  7.823848
675       <NA>        NA
676       <NA>        NA
677       <NA>        NA
678       <NA>        NA
679       <NA>        NA
680       <NA>        NA
681       <NA>        NA
682       <NA>        NA
683       <NA>        NA
684       <NA>        NA
685       <NA>        NA
686       <NA>        NA
687       <NA>        NA
688       <NA>        NA
689       <NA>        NA
690       <NA>        NA
691       <NA>        NA
692       <NA>        NA
693       <NA>        NA
694       <NA>        NA
695       <NA>        NA
696       <NA>        NA
697       <NA>        NA
698       <NA>        NA
699       <NA>        NA
700       <NA>        NA
701       <NA>        NA
702       <NA>  8.498721
703       <NA>  8.720931
704       <NA>  8.604290
705       <NA>  8.524848
706       <NA>  8.749027
707       <NA>  8.696090
708       <NA>  7.694210
709       <NA>  7.782043
710       <NA>  7.795972
711       <NA>  7.713504
712       <NA>  7.617341
713       <NA>  7.937261
714       <NA>  7.644573
715       <NA>  7.947938
716       <NA>  7.823768
717       <NA>  8.107298
718       <NA>  7.775563
719       <NA>  8.199186
720       <NA>  7.786336
721       <NA>  8.076256
722       <NA>  8.065975
723       <NA>  7.875243
724       <NA>  7.878596
725       <NA>  7.921734
726       <NA>  8.070790
727       <NA>  8.016949
728       <NA>  7.845682
729       <NA>  8.065453
730       <NA>  8.080804
731       <NA>  8.013844
732       <NA>  7.786325
733       <NA>  7.849526
734       <NA>  7.867606
735       <NA>  7.981242
736       <NA>  7.862458
737       <NA>  7.798821
738       <NA>  7.758436
739       <NA>  7.965872
740       <NA>  8.077637
741       <NA>  7.658674
742       <NA>  7.820142
743       <NA>  8.146252
744       <NA>  8.067008
745       <NA>  8.028474
746       <NA>  7.840641
747       <NA>  7.984290
748       <NA>  8.174601
749       <NA>  7.904545
750       <NA> 11.079144
751       <NA> 11.104733
752       <NA> 10.976559
753          Y 11.170351
754       <NA> 10.974355
755       <NA>        NA
756       <NA>        NA
757       <NA>        NA
758       <NA>        NA
759       <NA>        NA
760       <NA>        NA
761       <NA>        NA
762       <NA>        NA
763       <NA>        NA
764       <NA>        NA
765       <NA>        NA
766       <NA>        NA
767       <NA>        NA
768       <NA>        NA
769       <NA>        NA
770       <NA>        NA
771       <NA>        NA
772       <NA>        NA
773       <NA> 11.068500
774       <NA> 11.150093
775       <NA> 10.892873
776       <NA> 10.692017
777       <NA> 10.799624
778       <NA>  9.931291
779       <NA>  9.937099
780       <NA>  9.879677
781       <NA>  9.960879
782          Y  9.799230
783       <NA>  9.917411
784       <NA> 10.205820
785       <NA> 10.080766
786       <NA>  9.672058
787       <NA>  9.843718
788       <NA>  9.984209
789       <NA>  9.858062
790       <NA> 10.057493
791       <NA>  9.835000
792       <NA>  9.982838
793       <NA> 10.006017
794       <NA> 10.116749
795          Y  9.822760
796       <NA>  9.800464
797       <NA>  9.900778
798       <NA> 10.069179
799       <NA> 10.022430
800       <NA> 10.199161
801       <NA>  9.797960
802       <NA> 10.178221
803       <NA>  9.922937
804       <NA>  9.758456
805       <NA>  9.546896
806       <NA>  9.728604
807       <NA>  9.704398
808       <NA>  9.856627
809       <NA> 10.037302
810       <NA>  9.910142
811       <NA>  9.974996
812       <NA>  9.885505
813       <NA>  9.787729
814       <NA>  9.959380
815       <NA> 10.017679
816       <NA> 10.049657
817       <NA> 14.592119
818       <NA> 14.562083
819       <NA> 14.603564
820       <NA> 14.523114
821       <NA>        NA
822       <NA>        NA
823       <NA>        NA
824       <NA>        NA
825       <NA>        NA
826       <NA>        NA
827       <NA>        NA
828       <NA>        NA
829       <NA>        NA
830       <NA>        NA
831       <NA>        NA
832       <NA>        NA
833       <NA>        NA
834       <NA>        NA
835       <NA>        NA
836       <NA>        NA
837       <NA>        NA
838       <NA>        NA
839       <NA> 14.058246
840       <NA> 14.173888
841       <NA> 14.108883
842       <NA> 14.060153
843       <NA> 12.225855
844       <NA> 12.369737
845       <NA> 12.767654
846       <NA> 12.399036
847       <NA> 12.467906
848       <NA> 12.711704
849       <NA> 12.631612
850       <NA> 12.659693
851       <NA> 12.819407
852       <NA> 12.473868
853       <NA> 12.723541
854       <NA> 12.706043
855       <NA> 12.545120
856       <NA> 12.278360
857       <NA> 12.487621
858       <NA> 13.043044
859       <NA> 12.379962
860       <NA> 12.701159
861       <NA> 12.677753
862       <NA> 12.127693
863       <NA> 12.383393
864       <NA> 12.830199
865       <NA> 12.660316
866       <NA> 12.361482
867       <NA> 12.446998
868       <NA> 12.572500
869       <NA> 12.663130
870       <NA> 12.852704
871       <NA> 12.597568
872       <NA> 12.286835
873       <NA> 12.791770
874       <NA> 12.359842
875       <NA> 12.657701
876       <NA> 12.500388
877       <NA> 12.473674
878       <NA> 12.246329
879       <NA> 12.422503
880       <NA> 19.093727
881       <NA> 19.384591
882       <NA> 19.328623
883       <NA> 19.012358
884       <NA>        NA
885       <NA>        NA
886       <NA>        NA
887       <NA>        NA
888       <NA>        NA
889       <NA>        NA
890       <NA>        NA
891       <NA>        NA
892       <NA>        NA
893       <NA>        NA
894       <NA>        NA
895       <NA>        NA
896       <NA>        NA
897       <NA>        NA
898       <NA>        NA
899       <NA>        NA
900       <NA>        NA
901       <NA>        NA
902       <NA> 17.791539
903       <NA> 18.364952
904       <NA> 17.974789
905       <NA> 15.705892
906       <NA> 16.032357
907          N 15.868115
908       <NA> 15.962903
909       <NA> 15.718368
910       <NA> 16.028030
911       <NA> 15.788218
912       <NA> 16.146015
913       <NA> 16.017791
914       <NA> 16.094391
915       <NA> 16.059650
916       <NA> 15.974736
917       <NA> 16.089644
918       <NA> 15.904038
919       <NA> 15.730201
920       <NA> 15.945146
921       <NA> 16.019385
922       <NA> 16.017881
923       <NA> 15.728907
924       <NA> 15.809553
925       <NA> 15.936421
926       <NA> 16.064424
927       <NA> 16.131470
928       <NA> 16.031961
929       <NA> 15.992789
930       <NA> 16.071267
931       <NA> 16.103112
932       <NA> 15.955472
933       <NA> 15.986048
934       <NA> 15.763181
935       <NA> 15.696087
936       <NA> 15.896212
937          N 15.624137
938       <NA> 16.023563
939       <NA> 16.109588
940       <NA> 16.236257
941       <NA> 25.326674
942       <NA> 25.494612
943       <NA> 25.460209
944       <NA> 25.453316
945       <NA>        NA
946       <NA>        NA
947       <NA>        NA
948       <NA>        NA
949       <NA>        NA
950       <NA>        NA
951       <NA>        NA
952       <NA>        NA
953       <NA>        NA
954       <NA>        NA
955       <NA>        NA
956       <NA>        NA
957       <NA>        NA
958       <NA>        NA
959       <NA>        NA
960       <NA>        NA
961       <NA>        NA
962       <NA>        NA
963       <NA> 23.238457
964       <NA> 23.155656
965       <NA> 22.815140
966       <NA> 20.299463
967       <NA> 20.205599
968       <NA> 20.294681
969       <NA> 20.510868
970       <NA> 20.118232
971          Y 20.246755
972       <NA> 20.264198
973       <NA> 20.568773
974       <NA> 19.893345
975       <NA> 20.276619
976       <NA> 20.085955
977       <NA> 19.937756
978       <NA> 19.972539
979          Y 20.205227
980       <NA> 20.308143
981       <NA> 20.505487
982       <NA> 20.250556
983       <NA> 20.194189
984       <NA> 20.274195
985       <NA> 20.000511
986       <NA> 20.339655
987       <NA> 20.206581
988       <NA> 20.201431
989       <NA> 20.346296
990       <NA> 20.417110
991       <NA> 20.334573
992       <NA> 19.979562
993       <NA> 20.410593
994       <NA> 20.393707
995       <NA> 20.070019
996       <NA> 20.306643
997       <NA> 20.075175
998       <NA> 20.327450
999       <NA> 20.085347
1000      <NA> 20.468542
1001      <NA> 20.420857
1002      <NA> 33.167491
1003      <NA> 33.257160
1004      <NA> 33.088304
1005      <NA>        NA
1006      <NA>        NA
1007      <NA>        NA
1008      <NA>        NA
1009      <NA>        NA
1010      <NA>        NA
1011      <NA>        NA
1012      <NA>        NA
1013      <NA>        NA
1014      <NA>        NA
1015      <NA>        NA
1016      <NA>        NA
1017      <NA>        NA
1018      <NA>        NA
1019      <NA>        NA
1020      <NA>        NA
1021      <NA>        NA
1022      <NA>        NA
1023      <NA> 28.918686
1024      <NA> 25.192604
1025      <NA> 25.609985
1026      <NA> 25.343369
1027      <NA> 25.362742
1028      <NA> 25.378704
1029      <NA> 25.146247
1030      <NA> 25.274819
1031      <NA> 25.145572
1032      <NA> 25.574319
1033      <NA> 25.386432
1034      <NA> 25.638946
1035      <NA> 25.343586
1036      <NA> 25.654686
1037      <NA> 25.461572
1038      <NA> 25.160927
1039      <NA> 25.240468
1040      <NA> 25.725261
1041      <NA> 25.240549
1042      <NA> 25.376924
1043      <NA> 25.068693
1044      <NA> 25.505233
1045      <NA> 25.618353
1046      <NA> 25.454686
1047      <NA> 25.468733
1048      <NA> 25.147822
1049      <NA> 25.499008
1050      <NA> 25.416693
1051      <NA> 25.084569
1052      <NA> 25.474893
1053      <NA> 25.690995
1054      <NA> 25.230439
1055      <NA> 25.399971
1056      <NA> 25.331342
1057      <NA> 25.434679
1058      <NA> 25.233724
1059      <NA> 25.281693
1060      <NA> 42.406248
1061      <NA>        NA
1062      <NA>        NA
1063      <NA>        NA
1064      <NA>        NA
1065      <NA>        NA
1066      <NA>        NA
1067      <NA>        NA
1068      <NA>        NA
1069      <NA>        NA
1070      <NA>        NA
1071      <NA>        NA
1072      <NA>        NA
1073      <NA>        NA
1074      <NA>        NA
1075      <NA>        NA
1076      <NA>        NA
1077      <NA>        NA
1078      <NA>        NA
1079      <NA>        NA
1080      <NA>        NA
1081      <NA>        NA
1082      <NA>        NA
1083      <NA>        NA
1084      <NA>        NA
1085      <NA>        NA
1086      <NA>        NA
1087      <NA>        NA
1088      <NA>        NA
1089      <NA>        NA
1090      <NA>        NA
1091      <NA>        NA
1092      <NA>        NA
1093      <NA>        NA
1094      <NA>        NA
1095      <NA>        NA
1096      <NA>        NA
1097      <NA>        NA
1098      <NA>        NA
1099      <NA>        NA
1100      <NA>        NA
1101      <NA>        NA
1102      <NA>        NA
1103      <NA>        NA
1104      <NA>        NA
1105      <NA>        NA
1106      <NA>        NA
1107      <NA>        NA
1108      <NA>        NA
1109      <NA>        NA
1110      <NA>        NA
1111      <NA>        NA
1112      <NA>        NA
1113      <NA>        NA
1114      <NA>        NA
1115      <NA>        NA
1116      <NA>        NA
1117      <NA>        NA
1118      <NA>        NA
1119      <NA>        NA
1120      <NA>        NA
1121      <NA>        NA
1122      <NA>        NA
1123      <NA>        NA
1124      <NA>        NA
1125      <NA>        NA
1126      <NA>        NA
1127      <NA>        NA
1128      <NA>        NA
1129      <NA>        NA
1130      <NA>        NA
1131      <NA>        NA
1132      <NA>        NA
1133      <NA>        NA
1134      <NA>        NA
1135      <NA>        NA
1136      <NA>        NA
1137      <NA>        NA
1138      <NA>        NA
1139      <NA>        NA
1140      <NA>        NA
1141      <NA>        NA
1142      <NA>        NA
1143      <NA> 49.444588
1144      <NA> 49.831453
1145      <NA> 49.887663
1146      <NA> 50.063725
1147      <NA> 49.683286
1148      <NA> 49.748378
1149      <NA> 49.616993
1150      <NA> 49.675852
1151      <NA> 49.411498
1152      <NA> 50.062993
1153      <NA> 49.817167
1154      <NA> 49.895712
1155      <NA> 49.610599
1156      <NA> 49.719422
1157      <NA> 49.995702
1158      <NA> 49.759112
1159      <NA> 50.041707
1160      <NA> 49.958929
1161      <NA> 50.010466
1162      <NA> 49.619443
1163      <NA> 49.710237
1164      <NA> 49.859940
1165      <NA> 49.693797
1166      <NA> 49.666701
1167      <NA> 49.924587
1168      <NA> 49.932333
1169      <NA> 49.960272
1170      <NA> 49.809885
1171      <NA> 49.442550
1172      <NA> 49.986805
1173      <NA>        NA
1174      <NA>        NA
1175      <NA>        NA
1176      <NA>        NA
1177      <NA>        NA
1178      <NA>        NA
1179      <NA>        NA
1180      <NA>        NA
1181      <NA>        NA
1182      <NA>        NA
1183      <NA>        NA
1184      <NA>        NA
1185      <NA>        NA
1186      <NA>        NA
1187      <NA>        NA
1188      <NA>        NA
1189      <NA>        NA
1190      <NA>        NA
1191      <NA>        NA
1192      <NA>        NA
1193      <NA>        NA
1194      <NA>        NA
1195      <NA>        NA
1196      <NA>        NA
1197      <NA>        NA
1198      <NA>        NA
1199      <NA>        NA
1200      <NA>        NA
1201      <NA>        NA
1202      <NA>        NA
1203      <NA>        NA
1204      <NA>        NA
1205      <NA>        NA
1206      <NA>        NA
1207      <NA>        NA
1208      <NA>        NA
1209      <NA>        NA
1210      <NA>        NA
1211      <NA>        NA
1212      <NA>        NA
1213      <NA>        NA
1214      <NA>        NA
1215      <NA>        NA
1216      <NA>        NA
1217      <NA>        NA
1218      <NA>        NA
1219      <NA>        NA
1220      <NA>        NA
1221      <NA>        NA
1222      <NA>        NA
1223      <NA>        NA
1224      <NA>        NA
1225      <NA>        NA
1226      <NA>        NA
1227      <NA>        NA
1228      <NA>        NA
1229      <NA>        NA
1230      <NA>        NA
1231      <NA>        NA
1232      <NA>        NA
1233      <NA>        NA
1234      <NA>        NA
1235      <NA>        NA
1236      <NA>        NA
1237      <NA>        NA
1238      <NA>        NA
1239      <NA>        NA
1240      <NA>        NA
1241      <NA>        NA
1242      <NA>        NA
1243      <NA>        NA
1244      <NA>        NA
1245      <NA>        NA
1246      <NA>        NA
1247      <NA>        NA
1248      <NA>        NA
1249      <NA>        NA
1250      <NA>        NA
1251      <NA>        NA
1252      <NA>        NA
1253      <NA>        NA
1254      <NA>        NA
1255      <NA>        NA
1256      <NA> 63.469416
1257      <NA> 63.546102
1258      <NA> 63.523258
1259      <NA> 63.477376
1260      <NA> 63.961507
1261      <NA> 63.654355
1262      <NA> 63.628941
1263      <NA> 63.573938
1264      <NA> 63.880457
1265      <NA> 63.539050
1266      <NA> 63.540514
1267      <NA> 63.780502
1268      <NA> 63.856075
1269      <NA> 63.572275
1270      <NA> 64.000132
1271      <NA> 63.942538
1272      <NA> 63.625021
1273      <NA> 63.492469
1274      <NA> 63.707688
1275      <NA> 63.366794
1276      <NA> 63.471328
1277      <NA> 63.654416
1278      <NA> 63.573752
1279      <NA> 63.964729
1280      <NA> 63.752732
1281      <NA> 63.689620
1282      <NA> 64.071185
1283      <NA> 63.834859
1284      <NA>        NA
1285      <NA>        NA
1286      <NA>        NA
1287      <NA>        NA
1288      <NA>        NA
1289      <NA>        NA
1290      <NA>        NA
1291      <NA>        NA
1292      <NA>        NA
1293      <NA>        NA
1294      <NA>        NA
1295      <NA>        NA
1296      <NA>        NA
1297      <NA>        NA
1298      <NA>        NA
1299      <NA>        NA
1300      <NA>        NA
1301      <NA>        NA
1302      <NA>        NA
1303      <NA>        NA
1304      <NA> 64.774777
1305      <NA> 64.609002
1306      <NA> 64.875534
1307      <NA> 64.996694
1308      <NA> 64.750610
1309      <NA> 64.596887
1310      <NA> 64.849285
1311      <NA> 64.840499
1312      <NA> 64.707593
1313      <NA> 64.836489
1314      <NA> 64.688925
1315      <NA> 64.603284
1316      <NA> 65.038733
1317      <NA> 64.614364
1318      <NA> 64.710074
1319      <NA> 64.776279
1320      <NA> 64.592758
1321      <NA> 64.788182
1322      <NA> 64.702123
1323      <NA> 64.598391
1324      <NA> 64.351471
1325      <NA> 64.224045
1326      <NA>        NA
1327      <NA>        NA
1328      <NA>        NA
1329      <NA>        NA
1330      <NA>        NA
1331      <NA>        NA
1332      <NA>        NA
1333      <NA>        NA
1334      <NA>        NA
1335      <NA>        NA
1336      <NA>        NA
1337      <NA>        NA
1338      <NA>        NA
1339      <NA>        NA
1340      <NA>        NA
1341      <NA>        NA
1342      <NA>        NA
1343      <NA>        NA
1344      <NA>        NA
1345      <NA>        NA
1346      <NA> 65.315758
1347      <NA> 65.252923
1348      <NA> 65.271156
1349      <NA> 65.002114
1350      <NA> 65.373773
1351      <NA> 65.542051
1352      <NA> 65.069978
1353      <NA> 65.226503
1354      <NA> 65.096488
1355      <NA> 65.145541
1356      <NA> 65.397425
1357      <NA> 65.278940
1358      <NA> 65.422722
1359      <NA> 65.214001
1360      <NA> 65.196164
1361      <NA> 65.249922
1362      <NA> 65.492604
1363      <NA> 65.328235
1364      <NA> 65.179276
1365      <NA> 65.460647
1366      <NA>        NA
1367      <NA>        NA
1368      <NA>        NA
1369      <NA>        NA
1370      <NA>        NA
1371      <NA>        NA
1372      <NA>        NA
1373      <NA>        NA
1374      <NA>        NA
1375      <NA>        NA
1376      <NA>        NA
1377      <NA>        NA
1378      <NA>        NA
1379      <NA>        NA
1380      <NA>        NA
1381      <NA>        NA
1382      <NA>        NA
1383      <NA>        NA
1384      <NA>        NA
1385      <NA>        NA
1386      <NA> 65.111115
1387      <NA> 65.239517
1388      <NA> 65.491608
1389      <NA> 65.526845
1390      <NA> 65.750709
1391      <NA> 65.632116
1392      <NA> 65.521823
1393      <NA> 65.352583
1394      <NA> 65.040354
1395      <NA> 65.403154
1396      <NA> 65.424743
1397      <NA> 65.225789
1398      <NA> 65.147704
1399      <NA> 65.396397
1400      <NA> 65.275942
1401      <NA> 65.212436
1402      <NA> 65.317199
1403      <NA> 65.529183
1404      <NA> 65.358934
1405      <NA> 65.503301
1406      <NA>        NA

plot full predictions

  • What do you think if the predictions? Any better than model 1?
  • Do you agree with the Hmax values generated by the model?
  • If not, why do you think the model is generating those values? (also see next slide)
fit2.predictions %>% ggplot(aes(x=day)) +
  geom_line(aes(group=plantID, y=height_cm, color=pop), alpha=.2) +
  geom_line(aes(y=prediction, color=pop)) +
  scale_color_brewer(type="qual", palette = "Set2", guide="none") +
  facet_wrap(~pop)

DPR data + all predictions

  • This may be helpful in answering the last question
fit2.predictions %>% ggplot(aes(x=day)) +
  geom_line(aes(group=plantID, y=height_cm, color=pop),
            alpha=.2,
            data = {fit2.predictions %>% filter(pop=="DPR")}) +
  geom_line(aes(y=prediction, color=pop)) +
  scale_color_brewer(type="qual", palette = "Set2") 

Let’s try a new model

  • Having Hmax alone vary by population did not produce a very good fit.
  • Exercise: Create and fit a model that allows Hmax and k parameters to vary by population
  • Plot your results: Any better?
  • Perform a model comparison among all 3 models using the loo technique shown about: which fits best?

My formula

f3 <- brmsformula(height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta),
                  Hmin + delta ~ 1,
                  Hmax + k ~ pop,
                  nl=TRUE)

get the priors

Please review your notes and figure out how to get the possible priors for the formula on the previous slide

                prior class      coef group resp dpar nlpar lb ub       source
 student_t(3, 0, 4.2) sigma                                  0         default
               (flat)     b                           delta            default
               (flat)     b Intercept                 delta       (vectorized)
               (flat)     b                            Hmax            default
               (flat)     b Intercept                  Hmax       (vectorized)
               (flat)     b     popCC                  Hmax       (vectorized)
               (flat)     b    popDPR                  Hmax       (vectorized)
               (flat)     b                            Hmin            default
               (flat)     b Intercept                  Hmin       (vectorized)
               (flat)     b                               k            default
               (flat)     b Intercept                     k       (vectorized)
               (flat)     b     popCC                     k       (vectorized)
               (flat)     b    popDPR                     k       (vectorized)

Set the priors

What priors should we use?

prior3 <- c(set_prior("normal(5,3)", lb=0, nlpar="Hmin"),
            set_prior("normal(60,15)", nlpar="Hmax"),
            set_prior("normal(0,20)", nlpar="Hmax", coef="popCC"),
            set_prior("normal(0,20)", nlpar="Hmax", coef="popDPR"),
            set_prior("gamma(2,2)", nlpar = "k"), # gamma constrains to be > 0
            set_prior("normal(0, 1)", nlpar = "k", coef="popCC"),
            set_prior("normal(0, 1)", nlpar = "k", coef= "popDPR"),
            set_prior("gamma(20, 3)", nlpar="delta")

)

model summary

 Family: gaussian 
  Links: mu = identity; sigma = identity 
Formula: height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta) 
         Hmin ~ 1
         delta ~ 1
         Hmax ~ pop
         k ~ pop
   Data: datbcd (Number of observations: 907) 
  Draws: 4 chains, each with iter = 2000; warmup = 1000; thin = 1;
         total post-warmup draws = 4000

Population-Level Effects: 
                Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
Hmin_Intercept      4.04      0.52     3.01     5.05 1.00     3259     2279
delta_Intercept     6.67      0.53     5.66     7.76 1.00     2149     2221
Hmax_Intercept     64.96      1.38    62.33    67.78 1.00     2695     2120
Hmax_popCC        -14.57     20.87   -45.37    28.28 1.00     1680     2386
Hmax_popDPR       -28.01      8.58   -38.55    -5.96 1.00     2060     1122
k_Intercept         0.67      0.01     0.65     0.69 1.00     2363     2409
k_popCC             0.11      0.10    -0.03     0.34 1.00     1580     2591
k_popDPR            0.32      0.08     0.15     0.45 1.00     1987     1221

Family Specific Parameters: 
      Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sigma    11.17      0.27    10.67    11.72 1.00     4365     2709

Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
and Tail_ESS are effective sample size measures, and Rhat is the potential
scale reduction factor on split chains (at convergence, Rhat = 1).

model plot

pairs

Predictions

Model comparison

Next

  • delta + Hmax
  • delta + k + Hmax (doesn’t work well…)
  • Uncertainty on the predictions
  • On to random effects and all pops
  • CENSORED DATA?